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ABSTRACT 

We present interstellar absorption line measurements of the ions S III and Al III towards the stars 
P 1 Sco, \x Col, £ Per, £ Oph, p Leo, and HD 18100 using archival data from the Goddard High Resolution 
Spectrograph on board the Hubble Space Telescope. The ions Al III and S III trace heavily depleted and 
non-depleted elements, respectively, in ionized gas along the sightlines to these late-O/early-B stars. We 
use the photoionization equilibrium code CLOUDY to derive the ionization correction relating the ratio 
A(A1 III)/AT(S III) to the gas-phase abundance [Al/S] 4 (= log{A r (Al)/A^(S)}, -log{Al/S} Q ) in the ionized 
gas. For spectral types considered here, the corrections range from 0.1 to 0.3 dex and are independent 
of the assumed ionization parameter, i.e., the ratio of ionizing photon density to mean electron density. 

Using the results of these photoionization models, we find [Al/S]j ~ —1.0 in the ionized gas towards 
/3 1 Sco, £ Per, and C Oph; along the low-density path towards fi Col we find [Al/S]j ~ —0.8. Since S 
is not depleted onto grains these values of [Al/S],(~ [Al/H]j) imply that Al-bearing grains are present 
in the ionized nebulae around these stars. If the WIM of the Galaxy is photoionized by OB stars, the 
observations of p Leo and HD 18100 imply [Al/S].; = —0.4 to —0.5 in the WIM and thus the presence of 
dust grains containing Al in this important phase of the ISM. While photoionization appears to be the 
most likely origin of the ionization for Al III and S III, we cannot rule out confusion from the presence 
of hot, collisionally ionized gas along the sightlines to (3 1 Sco and HD 18100. We find that [Al/S], in 
the ionized gas along the six sightlines is anti-correlatcd with the electron density and average sightline 
neutral density. The degree of grain destruction in the ionized medium of the Galaxy is not much higher 
than in the warm neutral medium. The existence of grains in the ionized regions studied here has 
important implications for the thermal balance of these regions. 

Subject headings: ISM: abundances — dust, extinction — HII regions — Galaxy: halo — ultraviolet: ISM 



1. INTRODUCTION 

Warm (10 4 K) ionized hydrogen is an important com- 
ponent of our Galaxy's interstellar medium (ISM). The 
diffuse warm ionized medium (WIM) of the Galaxy has a 
mass surface density one third that of neutral hydrogen 
(H I), with an extended vertical scale height (h z « 900 
pc), and a power requirement equivalent to the total ki- 
netic energy injected into the ISM by supernovae (Kulkarni 
& Heiles 1987; Reynolds 1991b). The WIM has principally 
been studied through the detection of faint emission lines 
(e.g., Reynolds et al. 1998b) and through the dispersion 
measures of radio emission from distant pulsars (Taylor & 
Cordes 1993; Reynolds 1991a). 

The goal of this work is to determine if there is evidence 
for dust grains in the ionized medium of the Galaxy, in- 
cluding low-density H II regions and the diffuse WIM. The 
distinction between H II regions and the WIM in this work 
is mainly one of distance from the source of ionization and 
possibly fractional ionization. The existence of dust in 
the WIM of the Galaxy has important ramifications for 
the heating (and cooling) of the gas and possibly also for 
the power requirements of the ionization. Photoelectric 
emission of electrons from the surfaces of dust grains is an 
important source of heating in the warm neutral medium 
(WNM) of the Galaxy (Wolfire et al. 1995a). Reynolds & 
Cox (1992) have shown that, if grains are present in the 



WIM, photoelectric emission may be the dominant source 
of heating in the WIM and may also be responsible for 
the enhanced forbidden line strengths that are character- 
istic of the WIM (e.g., Dettmar & Shulz 1992; Rand 1997, 
1998). 

To provide evidence for the existence of dust in this 
important phase of the ISM, we will measure the gas- 
phase abundance of Al relative to S in ionized gas. In 
the warm neutral medium relative gas-phase abundances, 
derived from absorption line spectroscopy, have been used 
to infer the elemental composition of dust grains (see Sav- 
age & Sembach 1996b and references therein). We apply 
this method to the ionized medium of the Galaxy. 

Absorption line spectroscopy of interstellar material 
yields column densities of species independent of the pre- 
vailing physical conditions (e.g., T e , n e ). In Table [j] we 
list the properties of several possibly important probes 
of weakly ionized gas in the vacuum ultraviolet. Along 
with the wavelengths and /-values (from Morton 1991) 
for the transitions given, we also give relevant ionization 
potentials, logarithmic solar-system abundances of the el- 
emental species relative to H (Anders & Grevesse 1989; 
Grevesse & Noels 1993), and representative values of the 
gas-phase abundances of each of the elemental species in 
warm neutral halo gas (see Savage & Sembach 1996b). The 
lines and ionic species listed in Table |l| are those likely to 
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provide measurable column densities in the spectral range 
accessible to the UV spectrographs previously or currently 
on HST (~ 1120 - 3200 A) and lines accessible to the 
Far Ultraviolet Spectroscopic Explorer, which will cover 
the wavelength range ~ 905 — 1195 A after its launch 
in early- 1999. The last two columns of Table [l] give the 
expected column density of each ionic tracer for a com- 
pletely ionized region with halo-like abundances having 
log N (H II) = 19.0 (where N(H II) is expressed in cm~ 2 ) 
and total hydrogen density nn = 1.0 cm -3 , as well as the 
expected peak optical depth of each line in the case of pure 
thermal broadening with T = 10 4 K. These optical depths 
will be less if non-thermal broadening plays a significant 
role (a likely situation). 

Those ions from Table § accessible to HST include S III, 
Al III, P III, Si III, and Ti III. Given the great strength 
or weakness of many of the lines from Table l], one is typ- 
ically able to measure a tracer of a non-depleted element 
(see Spitzer & Fitzpatrick 1993; Fitzpatrick & Spitzer 
1994, 1997) in the ionized gas (S III) 2 and a tracer of 
a refractory element in ionized gas (Al III) in the HST 
bandpass. Relatively little is known about the Galactic 
distribution of these ionized gas tracers. Savage, Edgar, & 
Diplas (1990) used the International Ultraviolet Explorer 
to study the distribution of Al III in the Galaxy. They 
find an exponential scale- height of ft-Aini = 1-02^0 24 kpc 
for the distribution of Al III, thereby showing the distri- 
bution of Al III is slightly more extended than the neutral 
gas (hm = 0.67l^; 2 g kpc towards these same objects) and 
similar to that of free electrons. From the ratios of the 
estimated vertical column densities of Al III and free elec- 
trons (which are assumed to trace the protons of the WIM) 
they infer a gas-phase abundance of [Al/H]j > —1.7 in the 
ionized gas (using the solar system abundance Al/H given 
in Table||). Sembach & Savage (1992) have also studied 
the properties of Al III absorption towards a number of 
halo stars. They find that the absorption due to Al III is 
significantly different than that of the more highly-ionized 
species Si IV, C IV, and N V. 

In this work we present observations of interstellar ab- 
sorption due to the ions S III and Al III along the sightlines 
to six late-O/early-B stars: £ Persei (HD 24912), C Ophi- 
uchi (HD 149757), /3 1 Scorpii (HD 144217), \x Columbae 
(HD 38666), p Leonis (HD 91316), and HD 18100. These 
spectra are taken from the Goddard High Resolution Spec- 
trograph (GHRS) data archive. The stellar and sightline 
properties relevant to this study are given in Table Efl. The 
first four targets are relatively nearby disk stars (z <; 200 
pc). For these sightlines, the column density of ionized gas 
is likely dominated by H II region material, i.e., by fully- 
ionized gas immediately surrounding these stars. For the 
latter two stars, which are more distant and at high lati- 
tude, the paths through the ionized ISM are dominated by 
the WIM, by which we mean diffuse gas at large distances 
from any ionizing star. The species Al III and S III require 
more energy for creation than the H + typically used to 
probe the WIM; however, these ions likely represent sig- 
nificant fractions of the S and Al associated with the WIM. 
In all cases we depend upon the use of an ionization correc- 
tion factor to determine the relative abundances presented 



here. However, we will show that these small corrections 
are relatively insensitive to model assumptions. 

In §H we discuss the GHRS archival data and our re- 
ductions and analysis of these data. This includes in § [2~2] 
a description of the component-fitting analysis required 
for two of our six sigh tlin es, and measurements of Fe III 
from the literature in §2^3. We discuss the physical condi- 
tions and velocity structure of the ionized gas in §0 we in- 
clude in this a description of the complex kinematics along 
the sightline to £ Oph in §3.1.1. We use the photoioniza- 
tion code CLOUDY (Ferland 1996; Ferland et al. 1998) to 
model ionized regions about hot stars in §|] and use the 
results of these models to derive the gas-phase Al/S abun- 
dance in the ionized gas of the Galaxy. The abundances 
imply the existence of dust in t he io nized gas towards these 
six stars. We also discuss in §4.4 the confusing effects of 
collisionally ionized gas along our six sightlines. In §[5] we 
discuss the implications of our derived abundances and 
dust in the ionized phase of the Galaxy. A summary of 
our work and major conclusions is given in §|^. 

2. GHRS ARCHIVAL DATA 

In this section we discuss our reduction of the GHRS 
archival data and measurements of interstellar absorption 
features in these data. A more complete discussion can be 
found in How k, Savage, & Fabian (1998). Included in this 
section (§2.2) is a description of our approach to model 
component fitting. Component fitting is required to dis- 
entangle the S III A1190.208 A absorption line from the 
strong Si II line at Al 190.416 A in intermediate resolu- 
tion data. Also, we discuss the literature va lues of Fe III 
A1122.526 A absorption used in this paper (§|2.3|). 



2.1. Reductions and Measurements 

We have retrieved the GHRS archival data containing 
the S III and Al III transitions for the stars listed in Ta- 
ble |l|. The STScI archive identification codes and rele- 
vant information including exposure times, gratingmode, 
and aperture are given for each spectrum in Table @. The 
spectra have been calibrated using the standard CALHRS 
routine which includes conversion of raw counts to count 
rates and corrections for particle radiation contamination, 
dark counts, known diode non-uniformities, paired pulse 
events and scattered light. The final data reduction was 
performed using software developed and tested at the Uni- 
versity of Wisconsin-Madison. This includes the merging 
of individual spectra and allowing for additional refine- 
ments to the scattered light correction for the echelle-mode 
data (Cardelli, Ebbets, & Savage 1993). 

The spectra used for this study include data taken both 
before and after the installation of the Corrective Optics 
Space Telescope Axial Replacement (COSTAR). For in- 
formation regarding the pre-COSTAR performance of the 
GHRS, see Heap et al. (1995); the post-COSTAR per- 
formance of the GHRS is discussed by Robinson et al. 
(1998). The pre-COSTAR data used here were taken 
through the Small Science Aperture (SSA) , thereby avoid- 
ing the complication of the very broad wings found in the 
pre-COSTAR Large Science Aperture (LSA) line-spread 
function (LSF). The GHRS echelle-mode spectra (Ech-A 
and Ech-B) used here have a resolution of ^3.5 km s _1 



2 In many sightlines the SIR A1190.208 A line is blended with intermediate-velocity absorption from the strong Sill line at A1190.416 
(t) Ri +52 km s — 1 relative to SIR); this is particularly true for observations made at intermediate resolution (FWHM > 15 km s -1 ). 
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(FWHM). The observations of the S III A1190 A region 
using the G160M grating have a velocity resolution of ~ 20 
km s _1 , while in the Al III AA1855, 1862 A region the res- 
olution is ~ 11 km s" 1 (FWHM). For the G160M data we 
have used the SPYBAL observations taken before the sci- 
ence integrations to refine the absolute velocity scales (see 
Soderblom, Sherbert, & Hulbert 1993, 1994). The SPY- 
BAL observations contain emission lines observed in the 
wavelength range ~ 1500- 1540 A range. We have calcu- 
lated the average offset between the observed and expected 
wavelengths of these lines, and applied these wavelength 
offsets to the appropriate science observations (Soderblom 
et al. 1994). Using the SPYBAL observations, we find off- 
sets of —0.003 A and —0.033 A, respectively, are appropri- 
ate for the S III and Al III observations towards HD 18100. 
Towards (3 1 Sco the G160M observations of S III require an 
offset of —0.014 A. These wavelength adjustments should 
make the absolute velocities for the G160M observations 
good to ~ 3 km s . The velocity scale of the echelle data 
presented here is likely good to ±1 resolution element (i.e., 
3.5 km s- 1 ). 

Most of the data presented here employed the FP- 
SPLIT=4 procedure, which takes four exposures at differ- 
ent grating carousel positions. These scans are identified 
in Table ||. The FP-SPLIT exposures can be used to solve 
for the fixed-pattern noise spectrum of the GHRS (Cardelli 
& Ebbets 1994) and remove it from the observations. In 
general we have solved for the fixed-pattern noise, but only 
applied the solution in the cases where the noise spectrum 
showed significant structure within a few hundred km s _1 
of the line of interest. In this way we avoid adding noise to 
those spectra with few or no fixed-pattern features. Only 
for the ( Oph scans and the echelle observations of Al III 
towards /3 1 Sco was there structure deemed to be possibly 
significant to the current work; these are the only obser- 
vations for which we have applied the fixed-pattern noise 
solutions. 

We have normalized the spectra using low-order (< 4) 
Legendre polynomial fits to the local stellar continuum in 
regions free from interstellar absorption. Figures |l| and |^ 
present the normalized absorption line profiles of S III and 
Al III for the stars considered here. Also shown in Figure 
are the profiles of Zn II A2026.14 A or A2062.66 A. The 
Zn II profiles trace the column density of neutral material 
along these sightlines. Absorption from the Si II A1190 
A line can be seen in all of the plots showing the S III 
absorption. The expected shift of the Si II line relative to 
S III is +52.4 km s _1 . If intermediate negative velocity 
Si II is present, it can be blended with the S III line. The 
S III+Si II profile towards the star p Leo is a good ex- 
ample of this behavior; narrow intermediate velocity Si II 
features are present at vq w +12 and +25 km s^ 1 relative 
to the restframe S III along this sightlinc. These two fea- 
tures make it impossible to determine the S III absorption 
for vq > 5 km s -1 . Therefore our study of the p Leo sight- 
line is limited to vq < 5 km s . The behavior of Si II can 
be determined from the Si II A1193.29 A transition. 

Table || contains the measured equivalent widths and 
column densities of interstellar S III and Al III towards 
our six targets, along with the la uncertainties in these 
quantities. These uncertainties include contributions from 
photon statistics, continuum placement uncertainties, and 
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zero-level uncertainties. We have adopted a 2% zero-level 
uncertainty throughout. This may overestimate the errors 
in regions near heavily saturated lines, i.e., in the wave- 
length region of S III. See Sembach & Savage (1992; their 
Appendix) and Howk et al. (1998) for a more thorough 
discussion of these sources of uncertainty. 

For the echelle data the values of the column densities 
presented in Table ^ were derived by a straight integra- 
tion of the apparent column density profiles (Savage & 
Sembach 1991). In the absence of unresolved saturated 
structure, A a (A*) = N(X i ), where N a (X l ) and N(X i ) 
are the apparent and true column densities of the ion X' 1 , 
respectively. The S III lines observed towards £ Per and 
Q Oph are the only profiles for which we might expect the 
presence of unresolved saturated structure. Though the 
profiles seem resolved by the echelle-mode resolution, the 
value N(S III) for the £ Per and £ Oph sightlines may un- 
derestimate the total column density. For gas at T = 10 4 
K, the thermal Doppler spread parameters for S III and 
Al III are 6(S III) w 2.3 km s _1 and 6(A1 III) w 2.5 km s" 1 
(FWHM- 3.8 and 4.2 km s _1 , respectively). Thus the 
profiles of photoionized gas traced by these species should 
be marginally resolved, even without the likely addition of 
non-thermal broadening. 

As mentioned above, intermediate velocity Si II absorp- 
tion is seen towards p Leo at velocities near the restframe 
of S III. In both the Al III and Zn II profiles, an absorbing 
component is present near vq rs +20 km s _1 . The corre- 
sponding S III absorption is blended with the intermediate 
velocity Si II absorption. Therefore we report the column 
densities of Al III and S III for this sightline integrated over 
the velocity range Vq = —19 to +2 km s _1 . This only in- 
cludes the lower-velocity Al III absorbing component. We 
do not expect any Si II absorption in this velocity range 
given the profile of the Si II Al 193.29 A line present in the 
same GHRS observation. 

2.2. Component Fitting Measurements 

The S III Al 190.208 A lines observed towards /3 1 Sco and 
HD 18100 with the G160M grating are blended with the 
nearby Si II A1190.416 A line (v = +52.4 km s _1 relative 
to S III). To assess the extent of the blending, we have per- 
formed a component fitting analysis of these profiles using 
software kindly provided to us by E. Fitzpatrick and de- 
scribed in Spitzer & Fitzpatrick (1993; hereafter SF93). A 
model of the interstellar absorption was convolved with the 
instrumental LSF appropriate for the pre-COSTAR SSA, 
and the value of \ 2 minimized between this model and the 
data. The LSF adopted here for the SSA follows Spitzer 
& Fitzpatrick (1995) in using instrumental LSF given for 
the G160M grating in the GHRS Instrument Handbook 
(v3.0; Duncan 1992, Table 4-8). The LSF is well charac- 
terized by a Gaussian with a FWHM 1.08 diodes at the 
detector array. From these models, we derive the best-fit 
column density, Nk, Doppler parameter, bk, and central 
velocity, (v)k, for the k absorbing components of Si II and 
S III along these lines of sight. 

For the /3 1 Sco observations, we have used the the high- 
resolution GHRS Ech-B observations of the Si II A1808.013 
A line along this sightline to provide more information 
about the component structure of the Si II A1190 A line. 
We have fit the A1808 A line with a four-component model 
and subsequently fixed the values of Nk and bk from this 
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fit, though allowing common shifts in the velocities, to de- 
scribe the Si II A1190 A absorption. The S III column 
density presented for /3 1 Sco in Table |^ is derived from 
the best-fit one component model. Figure [l] displays the 
G160M data for this sightline as points, while the model 
derived from our fitting is displayed as a solid line. The 
central velocity for the S III profile is (vq) = —10.0 ± 0.4 
km s _1 (fitting error only). The 6- value is not well con- 
strained but is likely significantly less than the resolution 
of the G160M observations. 

For the sightline towards HD 18100, no high-resolution 
observations of Si II were available. However, we have fit 
the Si II A1193.290 A line present in the G160M obser- 
vation to assess the contribution of intermediate-velocity 
Si II absorption to the S III profile. The model fit to the 
Si II A1193 A profile is over-plotted as the dotted line on 
the S III profile for HD 18100 in Figure |. This model has 
been scaled by the appropriate factor to account for the 
differing oscillator strengths of the 1190 and 1193 transi- 
tions of Si II and shifted to the velocity scale of S III. It 
is clear from this presentation that there is a component 
of Si II overlapping the S III absorption in velocity. We 
have held the column density of this Si II component fixed 
when fitting the S III+Si II blend at 1190 A. We find the 
Si II component is relatively weak, of order 12% of the 
S III column density at these velocities. Table || gives the 
total column density derived for S III along this sightline. 
The total column density of S III we obtain agrees with 
the value log AT (S III) = 14.26 ± 0.06 derived by Savage & 
Sembach (1996a, hereafter SS96). 

We have similarly fit the Al III profile for HD 18100 
with a component model. We find evidence for at least 
two components in the S III and Al III absorption along 
this sightline. While the 6-values are ill-constrained and 
likely much less than the G160M instrumental resolution, 
the column densities and central velocities for the individ- 
ual absorbing components of S III are as follows: 

(u Q )i = -20±6 kms- 1 ; log N(S III) i = 13.51 ±0.25; 

(v Q } 2 = +5.8 ± 0.9 km s -1 ; log N(S III) 2 = 14.21 ± 0.03. 
For Al III we find: 

(u )i = -19±13 kms" 1 ; log JV(A1 IIl)i = 12.1 ± 0.4; 

(v G } 2 = +4.8 ± 1.6 km s" 1 ; log A(A1 III) 2 = 12.58 ± 0.11 

We do not believe unresolved saturation is a significant 
problem for either species for this sightline. Column densi- 
ties derived for the two Al III transitions through a straight 
integration of the are the N a (v) profiles are the same 
within the errors, and their N a (v) profiles show no evi- 
dence for saturation effects (see Savage & Sembach 1994; 
their Figure 4). For the S III profile we find that the equiv- 
alent width is within the range covered by the two Al III 
transitions, and the peak optical depth is very similar to 
that found in the Al III A1855 A profile (though it includes 
a small contribution from the overlapping Si II A1190 A 
line). The data for Al III and S III along the sightline 
to HD 18100 are shown in Figure || as points, while the 
component models for each species are plotted as solid 
lines. Considering the difficulty in deriving precise values 
for the (vq)i ~ —20 km s -1 component, we will use the 



integrated sight-line column densities in our analysis. We 
also note that the component structure buried in the inter- 
mediate resolution G160M data may be significantly more 
complex than our component-fitting analysis suggests (see 
SS96; Ryans, Sembach, & Keenan 1996). 

2.3. Literature Measurements of Felll 

For the stars /i Col, (3 1 Sco, and £ Oph measurements 
of the Fe III A1122.526 A line are available. Howk et al. 
(1998) report on the profile and column density of Fe III to- 
wards fi Col as measured with the GHRS using the G140M 
grating. Savage & Bohlin (1979) and Morton (1975) re- 
port on the gas-phase column density of Fe III towards 
(3 1 Sco and C, Oph, respectively, based upon Copernicus 
measurements. The column densities derived by these au- 
thors for Fe III are given in Table |3|. The data derived 
from the earlier Copernicus papers have been modified to 
reflect more recent determinations of the /-values of this 
transition (Morton 1991). We adopt / = 0.0788, which 
amounts to a modification of —0.15 dex to the column 
densities adopted in the works of Savage & Bohlin (1979) 
and Morton (1975). 

The Fe III absorption towards /3 1 Sco and £ Oph could 
be substantially blended with C I absorption at A1122.438 
A and A1122.518 A. Howk et al. (1998) have put restric- 
tive limits on the degree of contamination to the Fe III 
absorption towards \i Col. They find the level of contami- 
nation is not significant compared with the errors. For the 
Copernicus results, we have estimated the errors quoted in 
Table ||, though these column densities should be consid- 
ered upper limits due to the possible contamination from 
C I absorption. 

3. PHYSICAL CONDITIONS AND VELOCITY STRUCTURE 
OF THE IONIZED GAS 

The sightlines to the stars in our sample cover a wide 
range of physical conditions and environments. The ab- 
sorption lines due to ionized gas along the sightlines to- 
wards P 1 Sco, £ Oph, and £ Per are likely dominated by 
photoionized gas in the H II regions of these disk stars. 
The sightlines towards the higher-latitude stars sample the 
WIM of the Galaxy. Towards /i Col, the gas being sam- 
pled may well be a mixture of both H II region and WIM 
gas, though to our knowledge no H II region has previously 
been identified about this star. Given the range of condi- 
tions and velocity structure, we mention some of the more 
important aspects of these ionized gas properties here. 

3.1. H II Region Sightlines 

The H II regions of some of our disk stars have been 
studied in detail. Reynolds & Ogden (1982) and Reynolds 
(1988b) have studied the H II region S27 surrounding 
£ Oph. Their analysis of faint Ha and [S II] emission 
from the nebula suggests (n 2 ) 1 / 2 = Uef 1 ^ 2 « 4 cm -3 with 
T e w 6700 K, where n e is the local electron density and 
/ is the volume filling factor. Reynolds (1988b) has sim- 
ilarly studied Sivan 4, the H II region surrounding £ Per. 
The Ha and [S II] emission suggest (n 2 ) 1 / 2 = n e f 1//2 sa 1.4 
cm" 3 for this nebula with T e w 8000 K. Howk et al. (1998) 
and Shull & York (1977) have derived (n e ) = 0.2 cm -3 for 
the ionized gas along the line of sight towards fi Col, which 
is similar to that in the Galactic WIM. This average den- 
sity has been derived from analyses of the excited states 
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Sembach 1991) of the ions Al III, Si IV, and S III. The 
Si IV and Al III profiles have been scaled upwards by fac- 
tors of 40 and 100, respectively. The lower panel of Figure 
|H shows the log of the apparent column density ratios of 
Al III and Si IV to S III as a function of velocity. The 
value of log[iV(Al III)/AT(S HI)] varies by about one dex 
between the peak of the Al III profile at v© w —8 km s _1 
and the peak Si IV absorption near u© ~ —15 km s . 
The component at u© « — 15 km s _1 coincides in velocity 
with the dense neutral/molecular cloud seen in the Zn II 
profile (see also Savage, Cardelli, & Sofia 1992). 

Figure [| shows the complexities that arise when look- 
ing in detail at a given interstellar sightline. The other 
sightlines in our sample may be similarly complex, with 
complications hidden in t he absorption profiles. The mod- 
elling we present in §4.1 is highly idealized and does not 
account for the physics of stellar wind bow-shocks, for ex- 
ample. For the C Oph sightline we will give results derived 
from the integrated sightline values of logiV(AlIIl) and 
N(S III); we will, however, note the results for the case of 
log7V(AlIII)/7V(S III) w -2.0 along this sightline, which 
is appropriate for gas associated with the strongest com- 
ponent of Al III along this sightline. 

3.2. The High- Latitude Sightlines to HD 18100 and p Leo 

The average electron densities towards the higher- 
latitude stars can be estimated using the 2 Pz/2 fine struc- 
ture level of C II (denoted C II*). Savage & Sem- 
bach (SS96) have studied the excitation of C II* towards 
HD 18100. They find, assuming the excitation is caused 
by electron collisions, (n e ) — 0.071 cm~ 3 . This treatment 
assumes that all of the C II* and S II along the line of 
sight arise in the same gas and is a lower limit to the true 
electron density in the ionized gas. The value (n e ) for 
p Leo given in Table [l] is our own determination assuming 
electron collisions are populating the upper fine-structure 
level of C II. We have measured log N(C II*) = 14.13±0.02 
and logiV(SII) = 15.51 ± 0.02 from archival GHRS data 
for this sightline. The column of C II* was derived from 
observations using the G160M grating and should be con- 
sidered a lower limit given the possible presence of unre- 
solved saturated structure. The average electron density 
is then determined using Eqn. (7) of SF93, though we 
adopt [C/S] » [C/H] = -0.4 from Cardelli et al. (1996). 
We find (n e ) > (0.074 ±0.005) (T/6000K) 5 cm" 3 . Thus, 
(n e ) towards p Leo is similar to that towards the halo stars 
HD 18100 and HD 93521 as well as the direction towards 
3C 273, which also probes halo gas. 3 Reynolds (1991a) 
finds the high-latitude WIM is clumped in regions having 
electron densities (n e ) rj 0.08 cm~ 3 , a value quite close to 
the average densities derived for the high-latitude sight- 
lines discussed here. 

For the two distant high-latitude stars HD 18100 and 
p Leo we find a very good velocity correspondence between 
low-ionization gas (traced for example by Zn II) and ion- 
ized gas traced by Al III and S III. Savage & Sembach 
(1994, 1996) present GHRS observations of the high and 
low stages of ionization, respectively, towards HD 18100. 

3 The cooling rate per neutral H-atom for the p Leo sightline in the CII 158 fj,m line is therefore also similar to these other halo sightlines, 
with l c > (1.4 ± 0.3) X 10 — 26 ergs s — 1 H-atom — 1 . The cooling in the 158 /im line along the £ Oph and § Per sightlines may be as low as 
a factor of five below this value (Gry, Lequeux, & Boulanger 1992), further illustrating the differences in the ionized gas along the low- and 
high-latitude sightlines. 



of Si II and N II, respectively, and is a lower limit to the 
true electron density. 

The Ha emission towards £ Per is centered at vq = +7 
km s _1 (Reynolds 1988b). Emission from the nebula S220 
1° to the north, which is also thought to be powered by 
£ Per, is centered at i>© = +12 km s _1 . These velocities are 
well within the limits of Al III and S III absorption towards 
this star. Towards £ Oph the Ha emission is centered at 
vq = — 13 km s _1 (Reynolds 1988b), which is roughly con- 
sistent with the Al III absorption and the positive- velocity 
edge of the S III absorption. 

Recent observations with the Wisconsin Ha Mapper 
(WHAM) Fabry-Perot instrument (M. Haffner, private 
communication) have similarly shown that Ha-emitting 
gas towards p Col lies at velocities quite similar to the 
center of S III and Al III absorption along this sightline. 
The sightline towards this star may be associated with an 
H II region, the WIM along this sightline, or a combina- 
tion of both. At the distance and latitude of p Col, we 
expect only ~ 20% of the WIM in the direction of p Col 
to be in front of the star (Reynolds 1991b). Given that 
the Ha emission in this direction is at velocities consistent 
with the S III profile presented here, even though the WIM 
along this direction should mostly come from beyond the 
star, we will assume the emission and absorption along the 
sightline towards p Col are mostly probing H II region gas. 

When comparing the velocities of the tracers of ionized 
gas being studied here and tracers of neutral gas, such as 
Zn II, we see three distinct arrangements in the nearby 
disk stars. The ionized gas along the sightline to /3 1 Sco is 
found at velocities between those of the major components 
observed in the low-ionization species. The Al III and S III 
absorption towards p Col resides at velocities similar to the 
principal component seen in neutral gas tracers. For each 
of the stars ( Oph and £ Per we find a component of low- 
ionization absorption that is aligned with one edge of the 
S III profile. The low-ionization components found at this 
edge are significantly narrower than the S III absorption 
profile. There is no systematic trend when comparing the 
velocity structure of neutral and ionized gas tracers along 
these low-z sightlines. 

3.1.1. The Complicated Case of £ Oph 

The velocity structure observed in tracers of ionized gas 
along the path to C Oph is quite complex and warrants a 
more complete description. The profiles of Al III, S III, and 
Zn II towards C Oph shown in Figure [j] have distinctly dif- 
ferent velocity structure. The S III absorption is broader 
than either the Zn II or Al III profiles and is asymmetric. 
The Zn II absorption lines up well with the negative ve- 
locity edge of the S III profile, while the Al III absorption 
is aligned quite well with the positive velocity edge. 

Sembach, Savage, & Jenkins (1994) studied the Al III, 
Si IV, and C IV absorption along this sightline in some 
detail. These authors attributed the Al III and Si IV ab- 
sorption to the photoionized nebula about £ Oph and sug- 
gested the velocity offsets seen between these two species 
indicated the H II region was expanding. In Figure |^ we 
plot the apparent column density profiles (see Savage & 
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The low-ionization lines along this sightline are found at 
the same velocities as Al III and S III, as well as the high- 
ionization species such as Si IV and C IV; Savage & Sem- 
bach (1994) find the observed profile widths increase with 
the ionization potential of the species. 

Toward p Leo we also see that the Al III and Zn II 
profiles are quite similar, showing two principal absorbing 
components (blends) centered at velocities Vq rj — 7 and 
+20 km s _1 . To show this correspondence more clearly 
we plot the apparent column density profiles (Savage & 
Sembach 1991) of Zn II, S III, and Al III for the sightline 
towards p Leo in Figure ||. Also shown is the profile of the 
2 P 3 / 2 C II fine structure level as observed by the GHRS 
G160M grating with a resolution of ~ 17 km s" 1 (FWHM). 
The absolute velocity scale of the C II* observations was 
determined through the use of a SPYBAL observations as 
discussed in §|| The nominal wavelengths were shifted by 
-0.039 A (or -8.7 km s" 1 at A1335.77 A) based upon the 
analysis of the SPYBAL observations. 

The apparent column density profiles of Al III and S III 
have shapes very similar to the Zn II profile, and the align- 
ment of the Al III and Zn II component structure is excel- 
lent. The S III profile is systematically shifted by ~ —2.2 
km s _1 with respect to Zn II. This may be due to uncer- 
tainties in the absolute velocity scale; however, Fitzpatrick 
& Spitzer (1994) have found a similar offset (—2.4 km s _1 ) 
in their S III observations towards 7 2 Vel and suggested 
a possible error in the rest wavelength of the transition. 
Given the similarities of the two profiles, we have applied 
a +2.2 km s" 1 shift to the observed S III velocities in pro- 
ducing the profile in Figure [|. The profile for C II* , which 
is likely tracing thermal electrons that are responsible for 
exciting C II to the fine structure level (see SF93), also 
seems to trace the Zn II component structure very well, 
though the difference in resolution between the two data 
sets make a direct comparison somewhat uncertain. 

The detailed velocity correspondence between tracers of 
ionized and primarily neutral gas seen towards the two 
high-z stars in our sample suggests the ionized gas is as- 
sociated physically with the neutral material. There is no 
kinematic evidence that the ionized and neutral phases are 
spatially separated towards p Leo. The good velocity corre- 
spondence between neutral and ionized gas tracers could 
arise if S III and Al III absorption were tracing ionized 
edges of neutral clouds (e.g., McKee & Ostriker 1977) or if 
the clouds seen towards these high-latitude stars were par- 
tially ionized with neutral and ionized tracers mixed (e.g., 
SF93). The kinematic profiles in our data do not allow us 
to distinguish between these two scenarios, or other more 
complex arrangements of the two phases. The ionized and 
neutral media along the HD 18100 sightline seem to show 
a similarly close relationship, though we will not make as 
strong a claim in this case given the lower resolution of the 
data. The ionization produced by decaying neutrinos as 
proposed by Sciama (1995, 1997) would result in partially 
ionized gas. However, the decay photons from neutrinos 
with E 1 = 13.7 + 0.1 eV (Sciama 1995) are incapable of 
ionizing Al + and S + . Therefore other sources of ionization 
(e.g., star light, X-ray background photons, or collisional 
ionization) are required to produce the Al +2 and S +2 we 
observe. 

The close association of neutral and ionized tracers was 



also observed along the sightline to HD 93521 by SF93, 
a high-latitude star at z w 1500 pc from the Galactic 
plane. The component structure along the sightline to- 
wards HD 93521 is more complex than that towards p Leo, 
but the correspondence between C II* and tracers of neu- 
tral gas (e.g., S II) is very good. Spitzer & Fitzpatrick 
interpreted the data for the HD 93521 sightline to imply 
the free electrons and neutral gas were cospatial and well- 
mixed. Thus, they argue for the existence of a partially 
ionized phase of the ISM at high-z. 

Recent WHAM observations of the Perseus Arm have 
shown the intensity ratio of [O I] A6300 to Ha is in the 
range 0.01 to 0.04 (Reynolds et al. 1998a). The weakness 
of [O I] emission implies the fractional ionization in the 
WIM of the Galaxy is very high since the ionization of O 
and H are strongly coupled through charge-exchange reac- 
tions. Reynolds et al. comment that the observed ratio im- 
plies that most of the observed emission from the Galactic 
WIM cannot arise from partially-ionized gas along these 
sightlines, which probe gas at distances from the plane 
\z\ <; 300 pc. Similar results have been obtained for the 
WIM in the edge-on galaxy NGC 891 (Dettmar & Schulz 
1992; Rand 1998). The assumption by SF93 of a partially- 
ionized phase of the ISM would seem to be contradicted 
by these [O I] measurements, although the WHAM obser- 
vations were taken along a different path through the ISM 
and probe distances z <; 300 pc. 

Further high-resolution observations of neutral and ion- 
ized gas tracers will help to disentangle the connections 
between these phases of the ISM. Given the presence of 
C IV and N V absorption at similar velocities towards 
HD 18100, it is also possible that these high-latitude sight- 
lines are not tracing gas photoionized by stellar radiation, 
but rather photoionization by cooling hot gas or more com- 
plicated interactions between various phases of the ISM, 
such as conductive interfaces or turbulent mixing layers 
(see SS96 and references therein) . The contribution to the 
column densities of S III and Al III by hot, collisionally 
ioni zed gas could compromise our results and is discussed 



1.4 



Although the source of ionization of the WIM is not 
well understood, in §EO we will derive the gas-phase abun- 
dances in the ionized gas towards the high-latitude stars 
p Leo and HD 18100 by assuming the clouds to be pho- 
toionized by radiation from OB stars. 

4. GAS-PHASE ABUNDANCES IN IONIZED GAS 

Deriving the gas-phase abundance of Al relative to S 
from the column densities of Al III and S III requires the 
application of an ionization correction factor (ICF) to ac- 
count for the unobserved ionization stages of Al and S 
(primarily Al + and S + ). The logarithmic gas phase abun- 
dance normalized to solar, [Al/S]j, where the subscript i 
denotes this value in the ionized gas, is related to the mea- 
sured column densities of S III and Al III by 



[Al/S], ee log{V(Al +2 )/V(S+ 2 )} - log{Al/S} 
-log{x(Al+ 2 )/ a; (S+ 2 )}, 



(1) 



where a;(Al +2 ) ee Y(Al +2 )/V(Al) and x(S+ 2 ) = 
N(S +2 )/N(S) are the ionization fractions of A1+ 2 and S+ 2 
in the ionized gas. The ratio {x(Al +2 )/x(S +2 )}~ 1 is the 
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ICF, which we will write ICF(A1 +2 ). More generally in 
this work, ICF(X l ) = {x(X i )/x{S+ 2 )}- 1 . 

In this section we derive x(A\ +2 )/x(S +2 ) for moderate 
to low-density ionized gas near late-O/early-B stars. We 
then apply the ICFs to our measurements to derive gas- 
phase abundances [Al/S]j. 



4.1. CLOUDY Photoionization Equilibrium Models 

We use the photoionization equilibrium code CLOUDY 
(v90.04; Ferland 1996 and Ferland et al. 1998) to model 
the ionization and temperature structure of diffuse, low- 
excitation H II regions. Our models assume spherically 
symmetric nebulae excited by a single star. We use a 
volume filling factor /, which is the fraction of the vol- 
ume filled with constant hydrogen particle density nn (in 
cm~ 3 ); the rest of the space is assumed to be filled with 
very tenuous material and in these models is treated as a 
vacuum. From the radially-averaged ionization structure 
of the nebula, we can derive the ICF appropriate for di- 
rect measures of the column densities toward the central 
ionizing stars. 

We use ATLAS line-blanketed, LTE stellar atmosphere 
models (Kurucz 1991) as input spectra to the CLOUDY 
models. Our models follow the temperature and ioniza- 
tion structure of the model H II region from 0.3 pc dis- 
tance from the exciting star to the point where the elec- 
tron density falls to < 5% the ambient density riH (i-e., 
x(W) = 0.95). While the models we present here assume 
solar abundances with no dust opacity, we have found 
the inclusion of sub-solar abundances of the refractory el- 
ements and the addition of dust grain opacity does not 
significantly alter our conclusions. This is because highly 
refractory elements, i.e., Fe, Ni, Al, Cr, etc., are not dom- 
inant nebular coolants. Further, grain opacity tends to 
mimic the absorbing characteristics of H; including this 
opacity therefore removes the same ionizing photons as H 
(see Mathis 1986b). While the inclusion of photoelectric 
heating by dust and the absence of minor nebular coolants 
may have a more pronounced effect on the predicted emis- 
sion line strengths, these thermal differences cause little or 
no change in the predicted ionization fractions, and hence 
column densities, of the species we are considering. De- 
tailed descriptions of our use of CLOUDY models to gain 
information regarding the high-ionization species towards 
H Col and the contribution of H II region material to the 
neutral tracers towards this star can be found in Brandt 
et al. (1998) and Howk et al. (1998), respectively. 

We have computed a grid of model H II regions with 
varying input spectra and ambient densities. Table || 
gives the radially-averaged ionization fractions for sev- 
eral UV tracers of photoionized gas in model H II re- 
gions. These models used ATLAS model atmospheres with 
T e ff = 33, 000 K, logL*/ 'L@ = 4.4 and solar system abun- 
dances. This stellar effective temperature and luminosity 
are appropriate for fi Col or £ Oph. We have adopted 
/ = 1.0 in these models while varying the ambient density 
of the gas. Many authors define an ionization parameter, 
which is related to the ratio of hydrogen-ionizing photon 
density, n 7 , to particle density, tih, in the ionized region. 
A traditional definition of the ionization parameter (e.g., 



Mathis 1986b; Shields & Kennicutt 1995) is 

where L is the hydrogen ionizing photon luminosity 
in photons s^ 1 , and the Stromgren radius is Rs = 
[3 J L/(47rn2 [ /a B )] 1 / 3 , where cub is the case B recombina- 
tion coefficient of H (Ostcrbrock 1989), and / is again 

1/3 

the volume filling factor. Using this definition U cx n H , 
and we find 3U — (n 7 /nn), where the average is over vol- 
ume. In this work we will primarily adopt the equivalent 
definition given by Domgorgen & Mathis (1994; hereafter 
DM94) . We write the alternate ionization parameter q as 

Q = "h/ 2 ^50, (3) 

where L50 is the stellar ionizing luminosity in units of 10 50 
photons s _1 . This definition is related to the traditional U 
by q = 10 _50 (367rc 3 /a^)L T3 , though it removes the depen- 
dence on temperature that is hidden in as (Osterbrock 
1989). Bright high-density H II regions are typically de- 
scribed by models with log(g) ^> —1.0, while Domgorgen & 
Mathis use —4.0 <J log(g) <J —3.0 in modelling the Galac- 
tic DIG. We give the value of log(g) for each model in 
Table for comparison with the models of DM94, Mathis 
(1986a) and others. We also give the equivalent values of 
\og(U). The values U quoted assume as = 2.59 x 10~ 13 
cm 3 s _1 appropriate for T e = 10 4 K (Osterbrock 1989). 

In general the ionization and thermal structure of low- 
density models is determined by the ionization parameter, 
the adopted gas-phase abundances, and the shape of the 
ionizing continuum. Table ^ shows the ionization fraction 
of S +2 is a function of the ionization parameter q; how- 
ever the ratios of the ionization fractions x(Al +2 )/;r(S +2 ), 
a;(Si +2 )/x(S +2 ), and x(P +2 )/x(S +2 ) are quite insensitive 
to changes in q for a given stellar effective temperature. 
Further, these ratios, which give us the ICF for determin- 
ing the relative gas-phase abundances of Al, Si, and P to 
S, are not large. The greatest correction, for Al, is only 
of order - 0.2 dex. The ICF relating N(Fe III)/7V(S III) 
to [Fe/S]i is also relatively well-behaved with respect to 
the ionization parameter, but we shall see that it shows 
substantial variation with the input stellar effective tem- 
perature. 

Table I gives the ratios of ionization fractions for 
log(g) = —4.0 as a function of stellar effective temperature 
in the range 27,000 < T eff < 39,000. One can see that 
the relative ionization fractions x(X l )/x(S +2 ) are more 
sensitive to changes in the shape of the underlying contin- 
uum than to changes in the ionization parameter q. Even 
so, the spread of values x{A\ +2 )/x{S+ 2 ), x(Si +2 )/x(S+ 2 ), 
and a;(P +2 )/a;(S +2 ) are still relatively small. Table g also 
gives the average ionization corrections appropriate for 
relating the ratio of the various ionized species to their 
gas-phase abundances in the ionized gas. Two values 
are given: one representing the radially-averaged physical 
properties of the nebula (ICF/j ac ;) and one the volume- 
averaged properties (ICFy Q ;). The radially-averaged val- 
ues are appropriate for the sightlines that primarily probe 
absorption through H II regions towards the exciting stars. 
The volume-averaged ICFs are more appropriate for ran- 
dom interceptions of unrelated H II regions or the WIM. 
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The volume-averaged ICFs tend to weight the outer por- 
tions of the model nebulae more strongly, and therefore 
the volume-averaged values tend to favor lower stages of 
ionization. This can be seen in the difference between 
log(ICF(Si+ 3 ) v . o; ) and log(ICF(Si+ 3 ) i , ad ) in Table § the 
radially-averaged ICF suggests a greater fraction of Si is 
in Si +3 than in the volume-averaged case. The ICFs given 
at the bottom of Table || are averages over the seven stel- 
lar effective temperatures considered in our models with 
log(g) = —4.0 and —2.0. Also given are the standard de- 
viations about the means. The standard deviations give 
us a measure of how sensitive our derived ICFs are to 
uncertainties in the effective temperature of the ionizing 
radiation source, and we will use these values as estimates 
of the errors in our derived ICFs. 

In Figure || we show log[x(X l ) / x(S +2 )] as a function of 



T ef f for the ionized species X 1 = , Yc +Z , Si +2 , P+ 2 , 
C + , and N +2 . Also shown in the bottom panel of this fig- 
ure are the logarithms of x(S+ 2 ), x(A\ +2 ), and x(Fe +2 ). It 
is clear that the ICFs relating these ions to the gas-phase 
abundances of these elements (relative to S) in the ionized 
gas are functions of the assumed stellar effective tempera- 
tures. However, for Al III, Si III, and P III, the dependence 
on T e ff is relatively small, having a total spread of <J 0.25 
dex over the range 29,000 <T eff < 39,000 K. The bot- 
tom panel shows that even when S +2 and Al +2 are not 
the dominant ionization stages of Al and S, the fractional 
abundances of the two ions follow each other very well. 
The lcr uncertainties quoted for the ratios of ionic column 
densities using absorption line spectroscopy are often ^0.1 
dex (e.g., Savage, Cardelli, & Sofia 1992; SF93; Howk et 
al. 1998), which is similar to or greater than the standard 
deviations found for the predicted ICFs given in Table | 
for Al III, Si III, and P III relative to S III. 

The ICF for Fe III (also C III and N III) shows a greater 
dependence on the effective temperature of the underlying 
stellar atmosphere than those for Al III, Si III, and P III. 
However, even uncertainties on the order 0.2 to 0.3 dex for 
the value of ICF(Fe +2 ) can distinguish between depleted 
and non-depleted abundance ratios for an element that is 
typically found to be as heavily incorporated into grains 
as Fe. 

A concern when interpreting these results is the accu- 
racy of the atomic data used in deriving these models. 
In particular the atomic data for Fe are a concern (e.g., 
Pradhan & Bautista 1998). The atomic data adopted by 
CLOUDY are discussed in Ferland (1996) and Ferland et 
al. (1998). In general the photoionization cross-sections 
[using fits to the Opacity Project results by Verner et al. 
(1996)] and radiative recombination rates (see references in 
Ferland et al. 1998) are relatively secure. The dielectronic 
recombination coefficients ctdi for recombinations into the 
first two ions of Al are also relatively secure (Nussbaumer 
& Storey 1986). Though a dl (A\ +2 ) and a di (Al +3 ) are es- 
timated, Ferland (1996) predicts the ionization balance of 
Al° — Al +2 is relatively reliable. The low-temperature di- 
electronic recombination coefficients for the various ionic 
stages of S, Si, and Fe are not well constrained (Ferland 
1996; Ferland et al. 1998). 

CLOUDY estimates the unknown low-temperature di- 



electronic recombination coefficients of the first four ion- 
ization stages of elements in the third and fourth rows of 
the periodic table by adopting a mean of the rate coeffi- 
cients for the first four ionization stages of C, N, and O 
(Ali et al. 1991). 4 We have disabled CLOUDY's approxi- 
mation of ctdi for these elements to test its effects on the 
observed ratios. The ICFs derived for Al +2 , Si +2 , and P+ 2 
differ in models with and without the assumptions regard- 
ing dielectronic recombination by less than the standard 
deviations given in Table [s| For Fe III the models with 
ay^Fe 1 ) = 0.0 give values close to that given in Table 
|5[ However, the density dependence of a;(Fe +2 )/a;(S +2 ) 
is greatly increased. The values of ICF(Fe +2 ) is less cer- 
tain than the ICFs for the other species considered here 
given their strong dependence on the shape of the ioniz- 
ing spectrum. The same can be said for ICF(Si +3 ). The 
adopted error of a = 0.2 to 0.3 dex for ICF(Fe +2 ) likely 
encompasses enough of parameter space to make our error 
estimates not unreasonable. We will make estimates for 
[Fc/S]i where the data for Fe III are available, but it is 
important to recognize the strong dependence of the ion- 
ization correction on T e f f and the possibly inappropriate 
atomic parameters for Fe. 

4.2. [Al/SJi and [Fe/SJi in H II Regions 

In this section we derive [Al/S]i and [Fe/S]i for the sight- 
line s pr obing H II regions in the Galactic disk. The results 
of § [4.l| give ICFs that are directly applicable to H II re- 
gions about the stars £ Oph, £ Per, (3 1 Sco, and /j, Col. 
The distances to these stars are small enough that it may 
be reasonable to assume their H II regions dominate the 
column density of ions considered here. For £ Oph this 
comparison is complicated by the velocity structure, which 
suggests there are a number of processes at work along this 
sightline. A comparison of the integrated column densities 
of Al III and S III towards this star will not necessarily 
provide us a good meas ure of the gas phase abundance 
[Al/S],. As discussed in §3.1, the sightline towards \i Col 
may be probing H II region or WIM gas, or a mixture of 
both. We will assume here that the ionized gas along this 
sightline is probing an H II region. 

Table ^| presents the observed ratios of the integrated 
column densities of Al III and Fe III (where available) to 
S III for our program stars. Also given in this table are 
the derived logarithmic abundances [Al/S]i and [Fe/S]i. 
The abundances [Al/S]i and [Fe/S]i for the sightlines to 
[i Col, £ Per, (3 1 Sco, and C, Oph were derived using the 
ICFs appropriate for the stellar effective temperature s of 



these stars, with the error estimates as discussed in §4.1 



We are assuming in this approach that the measured ma- 
terial is mostly photoionized H II re gion gas. Given the 
velocity structure discussed in § |3.1.l| , [Al/S]; in the £ Oph 
H II region may be as high as —1.0. The values [Al/S]; 
and [Fe/S]i presented in Table ^ for these disk sightlines 
suggest grains are present in the H II regions surrounding 
all of these stars. 

Federman et al. (1993) present a GHRS measurement 
of P III A1334.813 A along the sightline to C Oph. They 
find logN(P III) = 12.94 ± 0.04. From the results pre- 
sented in Table |, we see logICF(P+ 2 ) = +0.01 ± 0.06. 



4 The adopted coefficients for the first four ionization stages of third and fourth row elements are: a di (X°) = 3 X 10 13 ; ayi (^ + ) = 3 X 10 
a<li(X +2 ) = 1.5 X 10~ u ; and a di (X+ 3 ) = 2.5 X 10~ n cm 3 s _1 . These are adopted from Nussbaumer & Storey (1983). 
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We thus find [P/S]; = -0.11 ± 0.07 for the ionized gas 
towards ( Oph. This value assumes the integrated sight- 
line column densities of P and S are tracing the same re- 
gions, which may be i ncorre ct given the velocity struc- 
ture towards this star (§3.1.1). The element P is thought 
to be very lightly depleted. In their study of the neu- 
tral gas abundances along the ( Oph sightline, Savage et 
al. (1992) find [P/H] = -0.23 for the warm neutral cloud 
centered at vq = —27 km s _1 (their component A). For 
the sightline towards p Col, Howk et al. (1998) derive 
[P/S] = — 0.03±0.03 in the low-velocity absorbing complex 
centered at v® = +23 km s _1 , which shows abundance 
patterns similar to component A towards £ Oph. The 
[P/S]j derived here for the C Oph sightline is consistent 
with previous [P/S] measurements for neutral material: P 
is lightly depleted if at all. 

4.3. [Al/S]i in the WIM 

In this subsection we consider the gas-phase abun- 
dances [Al/S]j in the halo WIM. The interpretation of the 
A^(A1 III)/A^(S III) measurements along the sightlines to 
p Leo and HD 18100 is subject to uncertainties regard- 
ing the ionization of the WI M and contamination from 
collisionally ionized gas (see §fl.4|). The long path- lengths 
through the WIM and the likelihood that any H II regions 
around these stars are in very low-density environments 
imply that a large fraction of the observed Al III and S III 
along these sightlines arises in the WIM of our Galaxy. 
Rho Leo should lie above ~ 60% of the WIM of the Galaxy, 
while HD 18100 should lie above almost all of the high- 
z layer of ionized material (Reynolds 1989; Savage et al. 
1990). Though the ionization source of the diffuse ionized 
material is not well known, the power requirements ioniza- 
tion by OB stars may be a viable mechanism (Reynolds 
1984; DM94). Here we estimate [Al/S]j towards p Leo 
and HD 18100 for the case where the WIM of the Galaxy 
is photoionized by OB stars. 

We will follow DM94 and model the WIM of the Milky 
Way using low-density H II region photoionization models. 
The ratio cc(Ai +2 )/ie(S +2 ) derived for our H II region mod- 
els is relatively insensitive to the ionization parameter and 
is not heavily dependent on the shape of the input ionizing 
spectrum; we therefor e be lieve it is appropriate to apply 
the results derived in §4.1 to the WIM of the Galaxy. Our 
observations of ionized gas towards p Leo and HD 18100 
probe the WIM and thus are not a simple radial inte- 
gration through an H II region. A volume-averaged set 
of ICFs are more appropriate for application to the WIM 
than the radially-averaged values used in the previous sub- 
section (D M94 ). Volume-averaged ICFs from the models 
treated in §4T are given in the last row of Table ||. We see 
that there are significant differences in the ICFs derived 
for most of the ions studied here when adopting volume- 
versus radially-averaged values. 

We do not know well the shape of the ionizing spectrum 
appropriate for the WIM. We will continue by adopting 
the mean of the volume-averag ed values ICF(A1+ 2 ) for 
our models. These values are given in Table | as ICF Vo! . 
The values [Al/S], given for HD 18100 and p Leo in Table 



|6| are therefore derived assuming log{a;(Al +2 )/x(S +2 )} = 
-0.37 ± 0.07, or logICF(Al+ 2 ) = +0.37 in the WIM. For 
the sightline towards p Col, the value [Al/S]j given in Ta- 
ble | was derived assuming logICF(Al +2 ) = +0.25 ± 0.07 
in an H II region about this star; if this sightline instead 
probes the WIM along this direction, the gas-phase abun- 
dance could be [Al/S^ = -0.66 ± 0.08. 

We can estimate values for logx(S +2 ) along our two 
high-latitude sightlines. Assuming N e — 7 x 10 19 / sin |6| 
cm~ 2 (Reynolds 1991b), and accounting for the 40% of 
this value expected to reside beyond p Leo, the predicted 
electron column densities along the high-latitude sightlines 
are log7V e = 19.7 and 19.9 for p Leo and HD 18100, re- 
spectively. Assuming N e w A(H + ) and the solar abun- 
dance of S/H (Anders & Grevesse 1989), for p Leo and 
HD 18100 we find logx(S+ 2 ) 1.3 and -0.9, respec- 
tively. Compared with our volume-averaged results for 
models having log(g) = —4.0, this implies the characteris- 
tic "effective temperature" for the ionizing radiation field 
of 28,000 £ T eff £ 33,000 K is appropriate for these 
sightlines. 5 This result is consistent with the limits derived 
from observations of diffuse He I recombination radiation 
(Reynolds & Tufte 1995). 

The Al III and S III data presented here, when com- 
pared with our model results, are inconsistent with a solar 
relative abundance of Al to S in the WIM. We have as- 
sumed that the intrinsic or cosmic abundance (gas+dust) 
of Al to S is given by the solar system values. Even in the 
presence of mildly sub-solar metallicities, [A/H] > —1.0, 
the abundance of Al relative to S is expected to be quite 
close to the solar abundance. Though Al is an odd-Z ele- 
ment, it is primarily produced in the C-,Ne-burning stages 
of massive stars. Thus, S and Al are both deposited into 
the ISM by Type II supernovae. 

In their study of elemental abundances in solar neigh- 
borhood low-mass stars, Edvardsson et al. (1993) derive 
abundances for Al as well as several a-elements. The be- 
havior of [Al/Fe] versus [Fe/H] in their dataset mimics 
that of the a-elements. 6 For their sample of 189 stars, 
all having [Fe/H] > —1.0, the average value of [Al/a] is 
+0.02 ± 0.07, where a in this context represents an aver- 
age over the elements Mg, Si, Ca, and Ti. The yield of 
Al relative to the a-elements does depend upon the ini- 
tial stellar abundances, but this effect is negligible for the 
range of metallicities of the Edvardsson et al. sample (and 
for our halo clouds). We do not expect the intrinsic ra- 
tio of Al/S, as long as the metallicity of the gas is ;> 0.1 
solar, to be significantly different than the solar-system 
value adopted here. Therefore the values of [Al/S]< < 0.0 
presented here for the WIM are not due to nucleosynthetic 
effects and imply the existence of dust in this phase of the 
ISM. 



4.4. Collisionally Ionized Gas and its Effects on [Al/SJi 

If hot, collisionally ionized gas is present along the sight- 
lines studied here, its imprint on the ionization balance of 
the gas could cause us to misinterpret the column density 

5 It is important to note that we are using the volume-averaged ICFs here, rather than the radially-averaged values given in Table pi For 
models having log (q) = -4.0 we find log x(S+ 2 ) = -1.38, -1.16, -0.99, -0.84, -0.77 and -0.53 for T eff = (27, 29, 31, 33, 35, 37, and 39) X 10 3 
K, respectively, in the volume- averaged model results. 

6 The element S is not included in their dataset, though see Wheeler, Sneden & Truran (1989) and references therein. 
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ratios given in Table |(| The collisional ionization equi- 
librium models of Sutherland & Dopita (1993), for exam- 
ple, suggest that for gas with temperatures T ~ 50,000 
K (log T = 4.7), log[a;(Al+ 2 )/x(S +2 )] = -0.93. There- 
fore, if one of our sightlines were dominated by gas in col- 
lisional ionization equilibrium at this temperature, with 
[Al/S]i = 0, i.e., no depletion, we would erroneously de- 
rive [Al/S]i « —0.5 by assuming the value of (ICF)voi 
given in Table |[ While the assumption of collisional ion- 
ization equilibrium is not likely to be valid, this calcula- 
tion suggests the presence of gas with T ^ 30, 000 K along 
our sightlines could compromise our conclusions. In this 
section we will give evidence suggesting, for most of our 
sightlines, it is unlikely our results are strongly affected by 
contamination from collisionally ionized gas. 

As a first step for approaching the problem, we have col- 
lected column density measurements of Si IV, and the ratio 
of Si IV and S III column densities, for all of our sightlines 
in Table jfl. Several of the Si IV column densities given 
in Table are taken from literature measurements using 
the GHRS (with references given in the table); others were 
derived from archival GHRS spectra using the techniques 
described in §2. With an ionization potential of 33.5 eV, 
the ion Si IV is potentially a tracer of collisionally ionized 
gas with T in the range (0.3 - 1) x 10 5 K, though Si IV 
can also be produced in warm gas via photoionization by 
hot stars. 

In the following discussion we consider the observational 
constraints regarding collisional ionization for each sight- 
line. 

C Op h - The absorption towards £ Oph, as discussed in 



|3.1.1 and shown in Figure shows evidence for changes 



in the ionization structure as a function of heliocentric 
velocity. The gas at vq <J —15 km s _1 is characterized 
by large amounts of Si IV absorption and relatively low 
amounts of Al III. Higher velocity gas along this sight- 
line, — 13<^vo<^— 5 km s _1 shows relatively little Si IV 
and higher ratios of Al III to S III. The peak of the Si IV 
profile may trace gas that is collisionally ionized. Our 
CLOUDY models, e.g., Table 0, cannot explain the large 
value of the ratio log[iV(Si JV)/N(S III)] « -1.8 observed 
near v Q w —20 km s^ 1 . However, Sembach et al. (1994) 
note that the Si IV line width implies gas with T < 5 x 10 4 
K. This temperature limit and the absence of associated 
C IV absorption lead Sembach et al. to suggest the Si IV 
absorption is due to photoionized gas in an expanding H II 
region. The gas between vq — —8 and —13 km s _1 is also 
likely photoionized gas in the H II region about £ Oph. 
This gas shows values of log[iV(Si YV)/N(S III)] <; -2.0, 
which can be generally explained by photoionization mod- 
els. The profile width of the Al III absorption measured 
by Sembach et al. (1994) implies gas with T < 3 x 10 4 
K. Morton (1975) also finds absorption from N II, N II*, 
N II** at i>0 = —8 km s~ 4 , suggesting the presence of 
(photo) ionized gas and a significant density of electrons 
at these velocities. Our earlier suggestion that a value of 
log[JV(Al III)/JV(S III)] « -2.0 is most appropriate for the 
photoionized H II region gas towards £ Oph is based upon 
the observed ratio at velocities corresponding to the peak 
of the Al III (and N II) absorption. Reynolds (1988b) has 
observed Ha emission centered at vq = — 13 ± 1 km s _1 
along this sightline, which is quite near the peak of the 
S III absorption. The temperature implied by an analy- 



sis of the breadth of the [S II] and Ha emission profiles 
is T = 6700 ± 700 K (Reynolds 1988b). Thus if we as- 
sociate the Ha and [S II] emission with the absorption 
seen in S III, Al III, and N II, the temperatures are too 
low to produce significant amounts of these ions through 
collisional ionization. 

£ Per - The sightline towards £ Per shows a com- 
plex of Ha-emitting regions at several different velocities. 
Reynolds (1988b) finds ionized regions excited by this star 
at velocities v@ = +3, +7, +12, and +14 km s _1 . Figure 
§ shows the N a (v) profiles of S III, Al III, and Si IV to- 
wards £ Per and the ratios of the latter two ions to S III, 
as Figure || did for £ Oph. The vertical dashed lines in this 
figure represent the velocities at which Reynolds (1988b) 
detects Ha emission, though shifted by —2.2 km s~ 4 . The 
peaks in the N a (v) profile for S III give a good match to 
the components detected in emission by Reynolds. In his 
analysis, Reynolds (1988) used the emission from [S II] and 
Ha to separate the thermal and non-thermal components 
of the observed velocity widths. He found that the emit- 
ting regions excited by £ Per had temperatures between 
T « 5,000 and 11,000 K, consistent with expectations 
for photoionized gas. We associate the emitting regions 
studied by Reynolds with the material seen in S III and 
Al III absorption in our GHRS spectra. At the implied 
temperatures, the contribution to the S III and Al III col- 
umn densities from collisionally ionized gas should be neg- 
ligible. The ICFs derived using photoionization models 
should therefore be appropriate for this material. 

/3 1 Sco - The derived log iV(Si IV)/N(S III) towards 
(3 1 Sco is quite similar to the measured values along the 
£ Oph and £ Per sightlines. Unfortunately, we have virtu- 
ally no information on the velocity structure of the S III 
absorption profile given the low resolution of the G160M 
data. The logarithmic ratio of log Y(Si IV)/iV(Al III) 
is in the range —0.2 to 0.0 over the velocity range of 
Al III absorption. The profile-weighted average veloci- 
ties (Sembach & Savage 1992) for Al III and Si IV are 
(vq) = —13.9 ± 0.6 and —11.4 ± 0.5 km s _1 , respectively. 
It is tempting to use the similarities between the C Oph 
and (3 1 Sco sightlines to argue that collisionally ionized gas 
is not confusing our analysis of this sightline. The neu- 
tral gas absorption profiles towards these stars are very 
similar (e.g., the Zn II profiles of Figure and the two 
stars lie at approximately the same distance from the sun 
and from the Galactic mid-plane. However, the differ- 
ences in spectral type would suggest that in the case of 
purely photoionized gas and equal gas densities, the ratio 
7V(Si IV)/iV(S III) should be smaller for (3 1 S co, i.e ., for 
the star with lower T e ff. Also, we have seen (§3.1.1) that 
there are some peculiarities along the ( Oph sightline, pos- 
sibly including the presence of collisionally ionized gas at 
velocities Vq <J — 15 km s~ 4 . We cannot rule out contribu- 
tions from collisionally ionized gas to the S III and Al III 
absorption towards /3 1 Sco. 

\x Col - For the sightline towards [i Col, measurements 
with the WHAM Fabry-Perot spectrometer (Reynolds et 
al. 1998) show the Ha emission is at velocities consistent 
with the S III absorption profile. The WHAM measure- 
ments, with a 1° beam, are well fit by a gaussian with 
(v Q ) = 22.9 ± 0.5 km s" 1 and b = 17.0 ± 1.5 km s" 1 
(M. Haffner, private communication). The breadth of the 
profile suggests that the gas traced by this emission has 
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T < 18, 000 K; in collisional ionization equilibrium models 
of gas at this temperature S III and Al III represent less 
than 1% of the total S and Al abundances (Sutherland & 
Dopita 1993). The profile- weighted average velocity (Sem- 
bach & Savage 1992) of the S III is (■%) = 23.2 ± 0.3 
km s _1 , with the Al III AA1855 and 1862 A lines having 
(v Q ) = 21.4 ± 1.0 and 23.9 ± 0.9 km s _1 , respectively. 
Due to the similarity of the velocities, we associate the 
Ha-emitting gas with the S III and Al III absorption. The 
temperature limits set by the WHAM spectra of this sight- 
line suggest that the contribution of hot, collisionally ion- 
ized gas to the S III and Al III column densities is small. 

p Leo - The archival GHRS spectra of Si IV towards 
p Leo show no evidence for absorption in the velocity 
range covered by the S III absorption studied here (i.e., 
-19 < v Q < +2 km s" 1 ). The 2a upper limit to iV(Si IV) 
given in Table @ implies log[7V(Si W)/N(S III)] < -2.1 in 
this velocity range, assuming &siiv ~ 10 km s _1 . Though 
there is Si IV absorption overlapping the higher-velocity 
Al III absorption, it appears to be significantly broader 
than the Al III absorption. We believe that the lack of de- 
tectable Si IV absorption at the velocities of the stronger 
absorbing complex seen in Al III (and S III), with limits 
that place the ratio N(Si IV)/AT(S III) at levels well below 
those detected along our H II region sightlines, suggests 
collisionally ionized gas is not an important contributor to 
the column densities of Al III and S III towards p Leo. 

HD 18100 - The sightline towards the star HD 18100 
has been studied by Sembach & Savage (1994), and their 
valu es of iV(Si IV) are given in Table [7]. As mentioned in 
§ |3.2| , the velocities of S III and Al III along this sightline 
are quite similar to those of the low ions (e.g., Zn II and 
Mn II), as well as those of the highly-ionized species Si IV 
and C IV. Indeed, Figure 4 of Sembach & Savage (1994) 
shows that the profiles of Al III, Si IV, and C IV are quite 
similar as viewed with the resolution of the GHRS G160M 
grating. Though the breadth of the profiles increases with 
ionization potential along this sequence, the increase seems 
mostly to occur towards negative velocities. It is possible 
that the highly ionized species discussed by Sembach & 
Savage have more complex velocity structures that would 
reveal themselves at higher resolution [the 1.4 km s _1 res- 
olution Ca II profile presented by SS96 and Ryans et al. 
(1996) shows evidence for many closely-spaced absorbing 
components along this sightline]. However, we have no 
evidence with which to rule out the possibility that the 
observed columns of S III and Al III are produced in col- 
lisionally ionized gas. The Si IV and C IV absorption 
profiles, which are at the same velocities as the lower ion- 
ization species, suggest that such contamination may in- 
deed be a real problem for this sightline. That the ratio 
of N(Si IV)/JV(S III) along this sightline is the highest in 
our sample also suggests collisionally ionized gas may be 
playing a role in determining the ionization balance in the 
gas we are studying. Thus our application of the ICFs 
derived from photoionization for this sightline may be in- 
appropriate. 

In summary it appears, with the exception of the sight- 
line to HD 18100 and possibly (3 1 Sco, collisional ionization 
does not likely explain the origin of the Al III and S III. In 
the case of HD 18100, the velocity overlap of the absorp- 
tion due to low, moderate, and highly ionized species, as 



well as the large column density of Si IV relative to S III, 
raises the possibility that much of the observed Al III and 
S III along this sightline could arise in collisionally ionized 
gas. 

5. DISCUSSION 

The Al III and S III column densities, when combined 
with the results of our photoionization modelling, suggest 
that Al-bearing dust is indeed present in the ionized gas 
along the sightlines considered here. For H II region gas 
this result is not surprising. The existence of dust in H II 
regions has been implied by several independent lines of 
reasoning (see Osterbrock 1989, Chap. 7; Mathis 1986b; 
Peimbert & Goldsmith 1972), though little can be dis- 
cerned about the composition of the grains or the degree 
to which dust grains are processed in the ionized ISM using 
these methods. 

Our measurements of [Al/S]j for the nearby sightlines 
that probe the H II regions surrounding the stars ( Oph, 
£ Per, (3 1 Sco, and p Col give us a measure of the in- 
corporation of Al into dust grains in these nebulae. As- 
suming [Al/S]j = [Al/HL, the dust-phase abundances 
of Al in the ionized material along these sightlines, in 
units of atoms per million hydrogen, are 10 6 (A1/H)rf = 
10 6 (Al/H) o -10 6 (Al/H), = 2.8 ±0.6, 2.6 ±0.5, 2.7 ±0.8, 
and 2.5 ± 0.5 for £ Oph, £ Per, (3 1 Sco, and p Col, respec- 
tively. These values are an order of magnitude less than 
the dust-phase abundance of Fe, Si, or Mg in the warm 
neutral ISM (Savage & Sembach 1996b; Howk et al. 1998). 
The dust-phase abundances of Al presented here are simi- 
lar to the values derived for Ni in the WNM, which has a 
similar solar system abundance (Savage & Sembach 1996b; 
Howk et al. 1998). 

Unfortunately, little is known about the gas-phase abun- 
dance of Al in the WNM of the Galaxy due to the great 
strength of the only available Al II transition at A1670 
A. Barker et al. (1984) have performed a curve-of-growth 
analysis of interstellar Al II absorption in the WNM us- 
ing data from the Copernicus observatory. However, for 
many of their sightlines the absorption due to Al is far 
up the flat portion of the curve of growth, making the 
uncertainties in the determination of the Al column den- 
sities quite large. In deriving the column density of Al II, 
Barker et al. employ an empirical curve of growth derived 
from Si II. However, since Al and Si are expected to have 
different gas-phase abundances from component to compo- 
nent due to changing depletion effects, it is likely that Si II 
and Al II have somewhat different curves of growth. Jenk- 
ins (1983), with several important caveats summarized at 
the beginning of his paper, has compared the equivalent 
widths of Al II A1670 A with the similarly strong Si II 
transition at A1304 A. While the nature of his comparison 
tends to heavily weight low-column density intermediate 
velocity features, Jenkins finds no evidence for changing 
[Al/Si] with z-height of the probe star with the implied 
value [Al/Si] w —0.2. This is near the upper range found 
by Barker et al. (1984), but is roughly consistent with their 
values. Again, there are many uncertainties with Jenkins' 
approach, among them that the Al II and Si II transitions 
are likely both on the flat part of the curve of growth, 
which may hide significant changes in the relative abun- 
dances. 

For comparison with our results, we look to those sight- 
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lines with the least saturated lines in the Barker et al. 
survey. For these sightlines Barker et al. consistently find 
[Al/H] w —1.1 to —1.0. This selection is strongly biased 
towards low-density sightlines. If we take [Al/S]j « —1.0 
as representative of the H II region sightlines in our sam- 
ple, this suggests the refractory grain material is not heav- 
ily affected by the conditions in low-density H II regions. 
However, the validity of this comparison is somewhat sus- 
pect given the uncertainties and biases in the Barker et al. 
results. 

Our estimates of [Fe/S]j in three of the sightlines con- 
sidered here (see Table yield values roughly consistent 
with those for [Fe/H]* in the Orion nebula from a variety 
of authors (Osterbrock et al. 1992; Peimbert et al. 1993; 
Baldwin et al. 1996; Rodriguez 1996; and Rubin et al. 
1997) and suggest Fe is incorporated into grains in these 
ionized regions. The relative values [Al/S]j versus [Fe/S]j 
suggest that the fraction of Al incorporated into grains is 
consistent with that of Fe or a bit less. For the sightlines 
towards ( Oph and fi Col we can directly compare the gas- 
phase abundances of Fe derived for the neutral and ionized 
sightlines. In each case, the velocity of the ionized gas is 
closest to the components that show the lowest value of 
[Fe/H], which are [Fe/H] = -2.4 for C Oph (Savage et al. 
1992) and [Fe/H] ss [Fe/S] = -1.3 for /i Col (Sofia et al. 
1993; Howk et al. 1998). Using the values of [Fe/S], given 
in Table [| we find that ~ 3% and ~ 9% of the dust-phase 
Fe in the neutral gas has been returned to the gas-phase 
in the ionized gas (with a factor of two uncertainty) . Our 
measurements of significantly sub-solar Al and Fe abun- 
dances in H II regions suggest the processing of grains in 
low-excitation H II regions is not much different than the 
destruction that occurs in the WNM of the Galaxy. 

The implications of dust within H II regions have 
been discussed by several authors (e.g., Sheilds & Kenni- 
cutt 1995; Henry 1993; McGaugh 1991; Aannestad 1989; 
Mathis 1986b). The most important effects are caused by 
the incorporation of possibly important coolants into the 
solid phase (Shields & Kennicutt 1995), the change in the 
thermal balance due to photoelectron emission (heating) 
and far- infrared thermal dust emission (cooling), and the 
competition of the dust opacity with H and He for ionizing 
photons (Mathis 1986b). McGaugh (1991) has shown that 
the existence of dust can significantly alter the Balmer line 
strengths expected from an H II region. This effect comes 
about because the dust is able to absorb a significant num- 
ber of photons that would otherwise go towards ionizing 
H (Mathis 1986b). As a consequence McGaugh suggests 
calculations of the star formation rates using only Balmer 
line intensities may underestimate the number of ionizing 
stars present. Also, the true volume of a dusty H II region 
will be smaller than the dust-free Stromgren volume. For 
low-density H II regions like those studied here, however, 
Mathis (1986b) has shown that the effects of dust absorp- 
tion may not provide a significant optical depth to ionizing 
photons. Mathis also points out that although dust is an 
additional source of opacity, the opacity due to dust does 
not affect the ionization balance of most species in H II 
regions given its similarity to the H opacity. Shields & 
Kennicutt (1995) discuss the important effects of dust on 
emission line strengths from metal rich [Z > Z®) H II re- 
gions, particularly those found near the centers of galaxies. 

Sembach & Savage (1996) have shown that in general, 



the gas-phase abundances of elements increases as one 
moves from the disk to the halo of our galaxy. This sug- 
gests an increasing degree of (incomplete) grain destruc- 
tion with increasing height above the plane of the Galaxy. 
In Figure || we show the gas-phase abundances [Al/S]j 
as a function of z-distance of the observed stars. Since 

5 is generally not depleted [Al/S], should closely follow 
[Al/H]j. There is a general trend of increasing gas-phase 
abundance of Al in the ionized gas with increasing height 
above the plane of the Galaxy. This is qualitatively con- 
sistent with the behavior observed by Sembach & Savage 
for the WNM. 

It is also known that the gas-phase abundances of el- 
ements increase with decreasing (uh) = A(H T)/d in the 
warm neutral ISM (Jenkins 1987; Savage & Bohlin 1979). 
Both the ionized and neutral gas densities are thought to 
decrease exponentially with z-height (Dickey & Lockman 
1990; Reynolds 1989), suggesting that the behavior seen in 
Figure^ may be tracing the density-dependence of the gas. 
In Figure ^ we plot the values [Al/S]j versus the electron 
density (top) and the average line of sight neutral density 
(bottom) , as given in Table H. Given the many definitions 
of the electron densities (average versus rms, etc.) we have 
used different symbols to represent the determinations of 
rms and average electron densities (see §||). 

Figure |t| shows a striking relationship between the av- 
erage densities and the gas-phase abundances of [Al/S]j 
along the sightlines considered here. We find that the gas- 
phase abundance [Al/S] j increases with decreasing electron 
densities. The observed relationship between [Al/S]j and 
n e is similar to the observed dependence of WNM abun- 
dances on average sightline neutral hydrogen (H I+H2) 
density (Edgar & Savage 1989; Jenkins 1987; Savage 

6 Bohlin 1979). The slope of [Al/S]* versus logn e is 
—0.37, similar to the value —0.38 derived for [Fe/H] versus 
log(nH) by Jenkins (1987). 

Though we are measuring the gas-phase abundance 
[Al/S]i, i.e., the abundance of Al to S in the ionized gas, 
the bottom panel of Figure M shows a significant correla- 
tion between [Al/S]* and logjn-a), the average line of sight 
neutral density. Savage et al. (1990) have also noted a cor- 
relation of log N(A\ III)/A(H I) with decreasing log(nn). 
These authors point out that this trend may be due to 
the changing ionization fraction of Al +2 with density, to 
changing values [Al/H]j with density, or perhaps both. 
Figure |?] shows that at least part of the trend observed 
by Savage et al. is due to the changes in the gas-phase 
abundance of Al with average neutral density. The trend 
observed by these authors, which is nicely matched by our 
data, and that seen in Figure ^ suggest that (nn) is a good 
indicator of the conditions in the ionized gas. The slope of 
[Al/S]j « [Al/H]j versus \og(nn) in our data is less steep 
than the slope of log A(A1 III) /N(R I) versus log(n H ) in 
the Savage et al. dataset, suggesting that a combination 
of changing gas-phase abundances and ionization fraction 
is causing the trend observed by Savage et al. and that 
perhaps the behavior seen in Figure |^ is more widespread 
than our six sightlines. 

In general Figure ^ implies that the ionized and neutral 
densities along a sightline follow the same trends, i.e., low 
(rin) also implies low values of n e . This relationship may 
simply be a manifestation of the known decrease in both 
ionized and neutral gas densities as a function of height 
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above the Galactic plane. This is also the expected be- 
havior if neutral clouds with ionized edges are providing 
most of the observed absorption (though see Reynolds et 
al. 1998a). Spitzer (1985) and Jenkins, Savage, & Spitzer 
(1986) have interpreted the dependence of elemental gas- 
phase abundances on the average sightline density (?ih) in 
the neutral ISM as a dependence on the relative contri- 
bution of two neutral media: clouds (cold and warm) and 
an intercloud medium. In this picture the warm intercloud 
medium has greater gas-phase abundances than the denser 
clouds. With an appropriate mix of each medium, the in- 
tegrated gas-phase abundance for a given line of sight can 
be reproduced. Perhaps a variation on this scenario is also 
appropriate for the ionized medium of our Galaxy. The 
relative mix of clouds may be similar between the neutral 
and ionized phases depending on the poorly-known rela- 
tionship between t hese phases. 

As discussed in §4.3, if the WIM of the Galaxy is pho- 
toionized by starlight from OB stars (DM94; Reynolds 
1984), our measurements of Al III and S III absorption 
towards HD 18100 and p Leo imply the existence of dust 
in this important phase of the ISM. The values [Al/S], for 
p Leo and HD 18100 are significantly higher than those 
for the other four stars (see Figures || and 0) , suggesting 
that the grain population in the high-latitude WIM has 
undergone a greater degree of processing (e.g., by shocks) 
than have the grains in the low-z H II regions. For the halo 
WIM we find a dust phase abundance of 10 6 (A1/H) d w 1.9 
to 2.1. It would appear that ~ 20% - 30% of the Al has 
been liberated from the solid phase in the high-z WIM 
compared with the results derived above for the disk H II 
regions. 

Very little is known about dust in the WIM of the Milky 
Way (or other galaxies) from previous studies. Though in 
principle one might be able to separate the thermal dust 
emission in the WIM from that of the H I and H2 gas, this 
has proven difficult. Boulanger et al. (1996) have studied 
the correlation of the far infrared (FIR) flux with the col- 
umn density of neutral hydrogen N(H I) at high Galactic 
latitudes. They find that the correlation between the A21- 
cm H I emission and the FIR emission detected by the 
Cosmic Background Explorer mission is quite good. They 
are not, however, able to rule out a dust abundance in the 
WIM similar to that observed in the neutral component. 

The existence of dust is important for maintaining the 
temperature of the WIM (Reynolds & Cox 1992; Dettmar 
& Schulz 1992). Reynolds & Cox (1992) show that the 
heating of the WIM may be in large part provided by pho- 
toelectron emission from grains. This requires the amount 
of grain heating per H atom to be similar to that found in 
the WNM. The presence of grains in the WIM has impor- 
tant ramifications for the diagnostic emission lines used to 
study this gas. Reynolds & Cox point out that the to- 
tal heating per H nucleus per second in the low-density 
WIM may be twice that of a typical, higher-density H II 
region. This has profound effects on the forbidden lines 
that provide the cooling for the gas. The increased cool- 
ing required over H II region gas increases the ratios of 
[S II], [N II], and [O III] to Ha over H II regions. Reynolds 
& Cox suggest this extra heating may in part be respon- 
sible for the enhanced forbidden line strengths observed 
from the WIM of our Galaxy and others (Reynolds 1985; 
Rand 1997). 



Our data suggest that the destruction of the dust grains 
in the high-z WIM has not been extreme. The value 
~ 20% — 30% given above for the amount of Al liberated 
from the solid- to gas-phase when going from disk H II re- 
gions to the halo WIM is consistent with differential mea- 
surements of warm neutral cloud abundances between the 
disk and the halo. For a small sample of halo and warm 
disk clouds, Howk et al. (1998) find roughly 20% - 30% of 
the dust-phase Fe and Si on average have been returned 
to the gas-phase between the disk and halo clouds. There- 
fore, the processing of grains in the halo WIM does not 
appear to be significantly greater than that experienced 
by clouds associated with the halo WNM. 

This discussion does not suggest that grain destruction 
mechanisms have not played an important role in the evo- 
lution of the gas being considered. The ~ 20% — 30% 
of dust-phase Al we see returned to the gas-phase in the 
halo WIM may be material that was initially bound in a 
refractory coating or mantle surrounding the grains. Sav- 
age & Sembach (1996b, see their Table 7) have tabulated 
the dust-phase abundances of a number of elements. The 
derived dust-phase abundances for the WNM of the Galac- 
tic disk give the composition of the grain cores and man- 
tles, while the observed abundances of the halo material 
gives information on the composition of the resilient grain 
cores that have probably been stripped of their mantles 
(Sembach & Savage 1996). Savage & Sembach (1996b) 
argue that the Fe returned to the gas-phase in halo ma- 
terial comes predominantly from the mantles thought to 
surround the resilient cores that survive the trip into the 
halo. They find the dust-phase abundance of Fe in grain 
cores to be 10 6 (Fe/H)c; = 25, while for the mantles they 
find 10 6 (Fe/H) rf = 7. Thus ~ 22% of the Fe incorporated 
into grains in the Milky Way disk resides in a mantle that 
is relatively easily stripped. The liberated Al seen in the 
WIM at high-z may also come from the mantles of grains, 
leaving the resilient grain cores to account for the remain- 
ing 70% — 80% of the Al missing from the gas phase. 

In considering the multiphase structure of neutral clouds 
in the Galactic halo, Wolfire et al. (1995b; see also Wolfire 
et al. 1995a) show that the stability of multiphase neutral 
clouds is affected by the intrinsic abundances in the gas 
and by the dust content of the clouds. Thus, the presence 
of dust in ionized halo clouds will have important impli- 
cations for the physical structure of the resulting neutral 
clouds if the ionized gas recombines. Indeed, in the case 
where the dusty multiphase clouds envisioned by Wolfire 
et al. (1995b) are situated above the disk of the galaxy, 
where they are bathed in ionizing radiation from the disk 
(assuming photons are able to leak out of the disk), they 
will be surrounded by ionized skins. If no processes beyond 
photoionization are responsible for producing the ionized 
edges of such clouds, the differences in the dust content of 
the neutral and ionized phases should be minimal. This 
may explain why the gas-phase abundances in the ion- 
ized gas towards p Leo and HD 18100, where the neutral 
and ionized phases of the ISM seem to coexist (at least in 
velocity space) , are so similar to the derived refractory el- 
ement gas-phase abundances in warm neutral halo clouds 
(Sembach & Savage 1996). 

6. SUMMARY 
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DUST IN THE IONIZED MEDIUM OF THE GALAXY 



This work represents one of the first absorption line 
studies of the WIM of the Galaxy. The observations im- 
ply the existence of Al- and Fe-bearing dust grains in the 
ionized gas of the Galactic disk and halo. 

A summary of the work presented here and our major 
conclusions is as follows: 

1) We present archival GHRS intermediate- and high- 
resolution absorption line observations of the moderately- 
ionized species Al III and S III in the ionized ISM towards 
six stars. The sightlines towards C Oph, £ Per, (3 1 Sco, and 
H Col probe primarily H II region gas in the Galactic disk. 
The extended high-latitude sightlines towards HD 18100 
and p Leo probe the WIM of the Galaxy at high z. Re- 
sults for Fe III from the literature are presented for p, Col, 
(3 1 Sco, and C Oph. 

2) We show, with the possible exceptions of the sight- 
lines to (3 1 Sco and HD 18100, that collisional ionization 
does not likely explain the origin of the observed amounts 
of Al III and S III. 

3) We have computed a grid of photoionization equi- 
librium models for low-density regions excited by late- 
O/early-B stars using the CLOUDY code (Ferland et al. 
1998). We show using our photoionization models that the 
ionization corrections for determining [Al/S];, [P/S]i, and 
[Si/S], using the species Al III, P III, Si III, and S III are 
relatively insensitive to the ionization parameter and the 
effective temperature of the ionizing spectrum. Deriving 
[Fe/S]i from the ratio N(Fe III)/7V(S III) requires a greater 
knowledge of the stellar effective temperature. 

4) We derive the logarithmic gas-phase abundances 
[Al/S]i « [Al/H]i in the ionized material towards the six 
stars in our sample using the results of our photoioniza- 
tion modelling. All of these stars have [Al/S]i ranging from 
— 1.2 to —0.4. Since S is normally not incorporated into 
dust, these abundance results indicate the incorporation 
of Al into dust grains in the ionized material along these 
six sightlines, though for the most distant stars, we cannot 
rule out the confusing effects of collisional ionization. For 
three stars we find [Fe/S]i ranges from —1.6 to —0.9. 

5) The gas-phase abundances [Al/S]i and [Fe/S]i de- 
rived here for the disk sightlines probing H II region ma- 
terial show that a significant degree of Al incorporation 
into grains is still present in the vicinity of stars, i.e., nei- 
ther the UV radiation fields from the stars nor any shocks 
associated with the stellar winds from these stars are suf- 
ficient to completely disrupt the refractory grains. The 
abundances we derive here are similar to the gas-phase 
Fe abundances [Fe/H]i derived for the Orion nebula using 
emission lines (e.g., Ostcrbrock et al. 1992; Baldwin et al. 
1996). 

6) If the WIM of the Galaxy is ionized by the light from 
OB stars (Reynolds 1984), the observed Al III and S III 
column densities towards HD 18100 and p Leo imply the 
existence of Al-bearing dust in the WIM. To our knowl- 
edge this is the first evidence for dust in the Galactic dif- 



fuse ionized gas, though the effects of collisional ionization 
cannot be ruled out, particularly for the sightline towards 
HD 18100. 

7) The gas-phase abundances [Al/S]^ in the ionized ma- 
terial increases with height z above the Galactic plane. 
Further, the values [Al/S]i increase with decreasing aver- 
age or rms electron densities and with decreasing average 
sightline neutral hydrogen density. This behavior is simi- 
lar to that of the gas-phase refractory abundances in the 
warm neutral medium of the galaxy (Jenkins 1987; Bohlin 
& Savage 1979). This general trend implies a greater re- 
turn of elements to the gas phase in more diffuse environ- 
ments. 

8) The observed values of [Al/S]i are similar to the abun- 
dances of other refractory elements seen in the WNM of 
the disk and halo. Further, the variation of [Al/S]i with 
density (electron or neutral) is also similar to the values 
observed for refractory elements in the WNM (esp., Fe 
or Mn). Our analysis implies that the processing of dust 
grains in the ionized gas may not be much different than 
that in the low-density warm neutral medium. 

9) We discuss the velocity structure of the WIM along 
the two high-latitude sightlines in our sample. These direc- 
tions show a very close relationship between the tracers of 
neutral material and ionized gas, similar to the correspon- 
dence observed by Spitzer & Fitzpatrick (1993). The data 
show no kinematic evidence for a separation of the ionized 
and neutral material. This is consistent with a partially- 
ionized medium in which the neutrals and ions are well 
mixed (e.g., SF93), the neutral clouds are surrounded by 
ionized envelopes (e.g., McKee & Ostriker 1977), or other 
more complex scenarios. 

We have benefitted from conversations with R. 
Reynolds, B. Benjamin, and J. Mathis in preparing this 
manuscript. We also thank E. Jenkins for a careful read- 
ing of our work and for suggestions that improved the final 
result. We thank E. Fitzpatrick for sharing his component 
fitting software with us, and M. Haffner, R. Reynolds, and 
the WHAM group for sharing their results in advance of 
publication. We feel G. Ferland and his co-workers at the 
University of Kentucky have done a great service for the 
astronomical community in making available their pho- 
toionization code CLOUDY. This work has only proceeded 
because of our confidence in the implementation of the 
physics included in CLOUDY. [The CLOUDY code and 
documentation, e.g., Ferland (1996), can be accessed v ia 



the address http : // www . pa . uky . edu/ ^gary/ cloudy/ 



the World Wide Web.J The archival GHRS data used 
here were acqui red from the GHRS Instrument De finition 
Team archive (http://hrssun.gsfc.nasa.gov/). This 
research has made use of the SIMBAD database, operated 
at CDS, Strasbourg, France. JCH recognizes support from 
a NASA Graduate Student Researcher Fellowship under 
grant number NGT-5-50121. 
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Table 1 

UV Probes of Weakly Ionized Gas 
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a Vacuum wavelengths from Morton (1991). 
b Oscillator strengths from Morton (1991). 

c The logarithmic "solar" abundances of the elements, log{A"/H} . We have adopted the solar system 
meteoritic abundances from Anders & Grevessee (1989) except for C which is the photospheric value from 
Grevesse & Noels (1993). 

d Typical values of the logarithmic normalized gas-phase abundance seen in warm neutral halo clouds, defined 
such that [A/H] = log{A(A)/A(H)} - log{A/H} . These values are for the warm cloud at v Q w +41 km s _1 
seen towards \i Col (Sofia, Cardelli, & Savage 1993; Shull & York 1977), which is typical of "halo"-type neutral 
clouds (Howk et al. 1998). 

e Expected column density of each species for a fully-ionized cloud with log A(HII) = 19.0 cm~ 2 as calculated 
using the CLOUDY photoionization equilibrium code (see §4.1). The assumed gas-phase abundances of each of 
the elements arc dictated by the solar system values modified by the normalized gas-phase abundances typical 
of halo material. The shape of the ionizing spectrum was taken to be that of an 09.5V star with T e ff — 33, 000 
K. 

f Expected peak optical depth of each line assuming the column densities given in column 8 with 6-values 
appropriate for gas at T = 10, 000 K and no non-thermal broadening. 

s The gas-phase abundance of C has only been measured accurately in cool neutral clouds found in the disk. 
We tentatively adopt the Cardelli et al. average [C/H] w —0.4 for our "halo" abundances. 

h The ratio Arl/HI has recently been measured to be significantly sub-solar along a number of low-iV(HI), 
partially ionized sightlincs by Sofia & Jenkins (1998). However, these authors argue that the large ionization 
cross-section of Ar° implies that much of the Ar may reside in the form of Ar + , which they did not observe. 
We therefore adopt [Ar/H] w 0.0 in warm neutral halo gas. 

1 The abundance of Al is poorly known in the warm neutral medium because the A1II A1670 A line requires 
very large saturation corrections (Barker et al. 1984). We have assumed [Al/H] « [Fe/H]. 
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Table 1 

Stellar and Sightline Properties 
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a Distances, with the exception of p Leo and HD 18100, are based upon Hipparcos measurements 
of stellar parallax (Perryman et al. 1997). The distance estimates to p Leo and HD 18100 are from 
Keenan, Brown, & Lennon (1986) and Diplas & Savage (1994), respectively. 

b Distance from the mid-plane of the galaxy given the derived distances. 

c Adopted effective temperatures for each of the stars. These data are from the following: Code et 
al. (1976; C Oph); Holmgren et al. (1997; /3 1 Sco); Howarth & Prinja (1989; p Col); Keenan et al. 
(1986; HD 18100); Keenan & Dufton (1983; p Leo); Sokolov (1995; 23 Ori); and Vacca, Garmany, & 
Shull (1996; £ Per). 

d The dust color excess, E(B-V), is taken from Diplas & Savage (1994). 

c Average line of sight neutral hydrogen densities towards these stars, where (n(H) } = {^V(H I) + 
27V(H 2 )}/d. We have taken values for iV(H I) from Diplas & Savage (1994) and N(R 2 ) from Bohlin, 
Savage, & Drake (1978). The value listed for HD 18100 is (n(H I)). 

f Electron densities along these lines of sight, when available. The data presented here are from Howk 
et al. (1998; p Col); Reynolds (1988b; £ Per and C Oph); Savage & Sembach (1996; HD 18100); and 
this work (p Leo). The values quoted for £ Per and £ Oph are the values (rig) 1 / 2 characteristic of the 
H II regions surrounding these stars. For p Col the value quoted is (n e ) of the H II region. The values 
quoted for HD 18100 and p Leo are the values (n e ) of the warm ionized gas in these directions. 
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Table 2 
Log of GHRS Archival Data 



Star 


Spectral Range 


Rootname a 


Exp. b 


Mode 


Aper. d 


FP-SPLIT/ 




[A] 




r i 

[secj 


& Order 




COSTAR? 6 


£ Per 


1 1 87 7-11 94 2 


Z0OY01 0LT 


691.2 


Frh-A /47 


SSA 


4/F 


fHD 9491 9) 

\ J_J_J__/ Li^-iJ J_ Li 1 


1857 1-18fi7 4 

lOJI .1 lOU 1 .4: 


Z0C1Y01 1 ST 

ZJUVJ J- UllkJl 


1 72 8 


FcVi-B /30 

J_J(__J.l J_J / JU 


SSA 


4 /F 


H Col 


1184.8-1191.1 


Z2AF010PT 


108.8 


Ech-A/47 


LSA 


0/T 


fHD 38666^ 

1 J.J.J--' UOUUU 1 


1 1 84 8-1 1 91 1 

J_ J.Ott.O J-J.C'J-. J- 


Z2C0090PP 

zJZj V^UUZ/U J- J- 


1 08 8 


F.ch- A /47 


T,SA 


IT 

U/ 1 




1846 7-1856 9 


Z2D401 18T 


54.4 


Erh-R /30 


LSA 


0/T 




1859.8-1869.8 


Z2CX010LT 


27.2 


Ech-B/30 


LSA 


0/T 




1855.9-1866.0 


Z2D4020KT 


54.4 


Ech-B/30 


LSA 


0/T 


COph 


1188.6-1195.0 


Z2VX010CT 


691.2 


Ech-A/47 


SSA 


4*/T 


(HD 149757) 


1189.8-1196.2 


Z2VX010ET 


691.2 


Ech-A/47 


SSA 


4*/T 




1856.9-1866.9 


Z0LD020TT 


172.8 


Ech-B/30 


SSA 


4*/F 


I3 l Sco 


1180.1-1216.3 


Z0YU010AT 


172.8 


G160M/01 


SSA 


4/F 


(HD 144217) 


1856.7-1866.7 


Z0YU020AT 


172.8 


Ech-B/30 


SSA 


4*/F 


p Leo 


1188.7-1195.3 


Z2ZX010CT 


870.4 


Ech-A/47 


SSA 


4/T 


(HD 91316) 


1852.7-1862.9 


Z0ZI0314T 


86.4 


Ech-B/30 


SSA 


0/F 




1853.5-1863.6 


Z0ZI0315T 


86.4 


Ech-B/30 


SSA 


0/F 




1854.3-1864.4 


Z0ZI0316T 


86.4 


Ech-B/30 


SSA 


0/F 




1855.0-1865.1 


Z0ZI0317T 


86.4 


Ech-B/30 


SSA 


0/F 


HD 18100 


1181.4-1217.6 


Z13Z010AT 


1324.8 


G160M/01 


SSA 


4/F 




1842.9-1876.9 


Z13Z010NM 


1209.6 


G160M/01 


SSA 


4/F 



a STScI archival rootname. 

b Total exposure time given in seconds. 

c Grating mode and spectral order used for the observation. 

d Aperture used for the observation. The LSA subtends 1'.'74 x l'/74 on the sky for post- 
COSTAR observations, 2'.'0 x 2'.'0 for pre-COSTAR data. The pre-COSTAR SSA subtends 0'/25 x 
Or 25, while the post-COSTAR SSA is 0722 x 0722 on the sky. 

e The number of FP-SPLIT sub-exposures composing each observation. Asterisks mark those 
observations for which we have explicitly derived the fixed-pattern noise spectrum and removed 
it. This column also notes with a "T" those observations taken after the installation of COSTAR, 
and with an "F" for those taken before COSTAR. 
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Table 3 

Equivalent Widths and Column Densities of S III, Al III and Fe III 







W x ± cr [mA] a 






logiV ± a [cm 2 ] 


b 


Star 


S III 
A1190.2 A 


Al III 

A1854.7 A 


Al III 

A1862.8 A 


S III 


Al III 


Fe III C 


£ Per 
It Col 

COph 

/3 1 Sco 
p Leo 
HD 18100 


87 ±3 
16.3 ±0.6 

66.6 ±1.4 
27±5 d 

12.7 ±0.8 
54±8 d 


15.0 ±2.0 

17.5 ± 1.6° 
74 ±6 


51.2 ± 1.7 
9.9 ± 1.5 

19.7 ±0.8 
9.5 ±0.7 

40 ±6 


14.82 ±0.02 
13.82 ±0.02 
14.76 ±0.02 
13.98 ±0.08 d 
13.72 ±0.03 
14.29 ±0.06 d 


12.85 ±0.02 
12.01 ±0.05 
12.42 ±0.02 
12.06 ±0.04 
12.06 ± 0.04 e 
12.70 ±0.04 


1 Q Q7 +0.09 

10.01 _o.ii 
13.45 ±0.10 
13.10 ±0.10 



a Equivalent widths W\ in mA for the lines of S III and Al III with la uncertainties. 

b Column densities of interstellar S III and Al III in units atoms cm~ 2 . Also given are the la error estimates for 
these measurements. 

c The Fe III column densities quoted here are taken from the following: Howk et al. (1998; p Col); Morton (1975; 
£ Oph); and Savage & Bohlin (1979; (3 1 Sco). These column densities are all derived from the Fe III A1122.5 A line. 
The latter two are based upon Copernicus observations, while the former comes from GHRS G160M observations. 
The Copernicus observations have been adjusted by —0.15 dex to account for newer oscillator strengths (Morton 
1991; / = 0.07884) and the errors are estimates by the current authors. In the case of £ Oph and possibly (3 1 Sco, 
C I absorption could be contributing to these column densities. For p Col Howk et al. (1998) have put restrictive 
limits on the degree of this contamination and find it not to be significant compared with the quoted uncertainties. 

d These values are based upon GHRS G160M data and have been derived through a component fitting analysis. 

°These values are for the velocity range —19 < v Q < +2 km s _1 , which corresponds to the uncontaminated 
velocity range for the S III absorption. The total integrated sightline values for Al III are W x = 29 ± 3 mA and 
logiV(Al III) = 12.27 ±0.05. 
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Table 4 

CLOUDY H II Region Model for an 09.5 V Star 



n H b log(g) c log(C/) d loga;(S +2 ) loga;(Y')/x(S +2 ) f 



[cm- 3 ] 








Al +2 


Fe +2 


Si+ 2 


Si+ 3 


p+2 


0.02 


-4.1 


-4.3 


-0.45 


-0.25 


0.24 


-0.15 


-2.03 


-0.01 


0.05 


-3.7 


-4.1 


-0.40 


-0.24 


0.21 


-0.13 


-2.08 


-0.01 


0.1 


-3.4 


-4.0 


-0.36 


-0.23 


0.19 


-0.12 


-2.14 


-0.02 


0.2 


-3.1 


-3.9 


-0.33 


-0.23 


0.17 


-0.11 


-2.22 


-0.02 


0.5 


-2.7 


-3.8 


-0.28 


-0.22 


0.14 


-0.10 


-2.36 


-0.02 


1.0 


-2.4 


-3.7 


-0.25 


-0.21 


0.12 


-0.09 


-2.46 


-0.02 


1.5 


-2.2 


-3.6 


-0.24 


-0.20 


0.11 


-0.09 


-2.55 


-0.02 


10.0 


-1.4 


-3.4 


-0.18 


-0.18 


0.08 


-0.07 


-2.94 


-0.02 


100.0 


-0.4 


-3.0 


-0.15 


-0.20 


0.10 


-0.06 


-3.77 


-0.03 



a The ionization fractions reported here are the radially-averaged values of nebular 
models with an input stellar effective temperature T e tf = 33,000 K and total lumi- 
nosity log L*/Lq — 4.4. The ionizing photon luminosity in these models is « 4 x 10 47 
photons s -1 . 

b Total hydrogen density (neutral plus ionized) used in the model. 

c q = nnf 2 L§o, where njj is the ambient hydrogen density, / the filling factor, and 
L50 the stellar ionizing flux in units of 10 50 photons s — 1 . 

d U = L/(ATTRgTiHc), where hh is the ambient hydrogen density and L is the 
hydrogen ionizing photon luminosity in photons s . The Stromgren radius is 
R s = [3L/(47m^fa B )] 1/3 , where aB is the case B recombination coefficient of H. 
We give values assuming as = 2.59 x 10 -13 cm 3 s _1 appropriate for T e = 10 4 K 
(Osterbrock 1989). 

c \ogx(X*)/x(S+ 2 ) = A/(S+ 2 )/A/(S T otai). 

f These columns give \ogx(X l )/x{S+ 2 ) = -loglCF (see §@). 
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Table 5 

CLOUDY H II Region Model for log(g) = -4.0 a 



[K] 


logx(S+ 2 ) 




log 


x(X*)/x(S+ 2 ) 








re 


Ol 


Ol 


T3 + 2 

r 


27,000 


-0.71 


-0.11 


0.55 


0.12 


-4.93 


0.12 


29,000 


-0.58 


-0.16 


0.40 


-0.01 


-3.85 


0.05 


31,000 


-0.49 


-0.21 


0.30 


-0.10 


-2.78 


0.01 


33,000 


-0.44 


-0.25 


0.23 


-0.14 


-2.04 


-0.01 


35,000 


-0.40 


-0.30 


0.15 


-0.19 


-1.58 


-0.04 


37,000 


-0.37 


-0.33 


0.10 


-0.22 


-1.43 


-0.06 


39,000 


-0.34 


-0.36 


0.04 


-0.27 


-1.30 


-0.10 


log 


(ICF) flad b 


+0.24 ±0.07 


-0.17 ±0.19 


+0.10 + 0.10 


+2.6 + 1.4 C 


+0.01 + 0.06 


log(ICF) yoi d 


+0.37 ±0.07 


-0.4 ±0.3 


+0.38 + 0.21 


+4.3 + 1.4 C 


-0.05 + 0.08 



a The ionization fractions and ratios reported here are the radially-averaged values of nebular 
models with log(g) = log(nH,/ 2 i5o) = —4.0. The alternate definition of the ionization parameter 
gives \og{U) = log{i/(47ri?|nHc)} = —4.22 assuming a.B — 2.59 x 10~ 13 cm 3 s" 1 appropriate for 
T e = 10 4 K (Osterbrock 1989). 

b The mean value of the radially-averaged ionization corrections and standard deviation about 
the mean for each species considered. These numbers are the average and dispersion of the ICFs 
for all stellar effective temperatures considered for models having log(q) = —4.0 and —2.0. 

c The values of ICF(Si +3 ) are highly uncertain given the strong dependence upon the stellar 
effective temperature and ionization parameter. 

d The mean value of the volume-averaged ionization corrections and standard deviation about 
the mean for each species considered. These numbers are the average and dispersion of the ICFs 
for all stellar temperatures considered for models having log(g) = —4.0 and —2.0. 



Table 6 

Column Density Ratios and Logarithmic Gas-Phase Abundances 



Star 


logJV(Al m)/N(S III) 


logJV(PeIII)/JV(S III) 


[Al/S]i 


[Fe/S], 


£ Per 


-1.97 + 0.03 




-0.87 + 0.08 




fj, Col 


-1.82 + 0.05 


-0.40 + 0.11 


-0.78 + 0.08 


-0.87 + 0.21 


C Oph a 


-2.23 + 0.03 


-1.05 + 0.10 


-1.19 + 0.08 


-1.52 + 0.21 


(3 1 Sco 


-1.92 + 0.09 


-0.88 + 0.12 


-1.00 + 0.11 


-1.59 + 0.21 


p Leo 


-1.66 + 0.06 




-0.50±0.10 b 




HD 18100 


-1.59 + 0.07 




-0.43±0.10 b 





a The values given here for £ Oph are the for the integrated sightlinc. 

b The values given here for [Al/S]i towards p Leo and HD 18100 assume the shape of the 
ionizing spectrum of the diffuse ionized gas is not very different than that of a star with 
27, 000 ^i;^ 39, 000 K. 
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Table 7 

Measurements of Si IV Column Densities 



Star 


logA(Si IV) 


log[A(Si IV)/JV(S III)] 


Ref. 


£ Per 


12.89 ±0.03 


-1.93 


1 


p Col 


12.17 ±0.05 


-1.65 


2 


COph 


12.79 ±0.02 


-1.97 a 


3 


I3 1 Sco 


12.02 ±0.03 


-1.96 


1 


p Leo 


< 11. 6 b 


< -2.1 b 


1 


HD 18100 


13.10 ±0.04 


-1.19 


4 



References.— (1) This work; (2) Brandt et al (1998); (3) 
Sembach et al (1994); (4) Sembach & Savage (1994). 

a The value logLV(Si IV)/iV(S III)] given for the C Oph sight- 
line is the value integrated over all velocities. Figure || shows 
that this ratio is a strong function of velocity. Near v Q ss —15 
km s _1 , corresponding to the peak of the Si IV absorp- 
tion, a value of —1.6 to —1.8 is more appropriate; while at 
the peak of the Al III absorption near vq « —8 km s , 
logLV(Si IV)/iV(S III)] w -2.4 is more appropriate. 

b This 2cr upper limit to log iV(Si IV) for the p Leo sightlinc 
only applies to the velocity range vq = — 19 to +2 km s _1 (the 
range over which the S III measurements are made). There is 
detectable Si IV absorption at more positive velocities. The in- 
tegrated sightline column density is log AT(Si IV) = 12.24+o.og. 
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m o ^ — i lo o lo o 1 



Fig. 1. — Continuum-normalized absorption line profiles of Zn II, Al III, and S III for the disk stars £ Per, ( Oph, and /3 Sco are displayed 
on a heliocentric velocity scale. Echelle-mode data are plotted as histograms. The G160M data for the S III profile towards /ii_Sco are plotted 
as points, while the component model convolved with the instrumental spread function is overplotted as the solid line (sec {2.2). Strong Si II 
A1190.416 A absorption is seen in the S III A1190.208 A region of the spectrum (v = +52.4 km s _1 relative to S III). 
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Fig. 2.— As for Figure |j but for the higher-z stars \i Col, p Leo, and HD 18100. The G160M data for the Al III and S III profiles towards 
HD 18100 are overplotted with the component model convolved with the instrumental LSF; thti-dashed line in the S III profile for this star 
shows the absorption model for Si II A1190 A derived from the Si II A1193 A transition (see §2.2). 
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Fig. 3. — The top panel shows the apparent column density profiles of Al III, S III, and Si IV observed towards £ Oph at 3.5 km s — 1 
resolution displayed on a heliocentric velocity scale. The profiles of Al III and Si IV have been scaled upwards by factors of 100 and 40, 
respectively. The bottom panel shows the logarithm of the ratio of Al III and Si IV to S III as a function of velocity. 
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Fig. 4. — Apparent column density profiles of S III, Al III, C II*, and Zn II towards p Leo arc plotted against heliocentric velocity. The 
C II* observations were obtained with the Gf60M grating, the other measurements were made with the GHRS echelle gratings. The profiles 
of the first three ions have been scaled to match that of Zn II. The scale factors are noted in the plots. The S III profile has been shifted 
by +2.2 km s^ 1 (see text). The S III profile is not plotted for vq > 5 km s _1 because of strong contamination with intermediate negative 
velocity Si II Alf90.4f6 A absorption components. 
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Fig. 5. — The values of log x(X t ) / x(S +2 ) = — logICF(X l ) for several tracers of ionized gas as a function of the effective temperature, T e ff, 
of the ionizing star. The fractional ionization x(X % ) = N(X Z )/N(X) is derived using column densities integrated from the star to the edge of 
the model H II region. These data are for models characterized by the ionization parameter log(<jr) = —4.0. Also shown in the bottom panel 
is the behavior of the ionization fractions i(S+ 2 ), x(Fe+ 2 ), and x(Al+ 2 ) as a function of the assumed stellar effective temperature. 
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Fig. 6. — The behavior of [Al/S]i = log(Al/S)i — log(Al/S)Q, as a function of z-height of the probe star. Since S is not depleted onto grains 
[Al/S], ~ [Al/H]j gives a measure of the fraction of Al incorporated into dust grains in the ionized gas. 
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Fig. 7. — The behavior of [Al/S]; = log(Al/S); — log(Al/S)0 as a function of average density. Since S is not depleted onto grains 
[Al/S], ~ [Al/H]j gives a measure of the fraction of Al incorporated into dust grains in the ionized gas. The bottom panel shows [Al/S], 
versus the average sightlinc neutral density (iih) = [iV(H I) + 27V(H2)]/(i The top panel shows [Al/S]j versus the electron density n e . The 
filled circles in the top panel represent determinations of rms electron densities {n^} 1 / 2 , while the triangles are for determinations of average 
electron densities (n e ). 
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Fig. 8. — The top panel shows the apparent column density profiles of Al III, S III, and Si IV observed towards £ Per at 3.5 km s _1 
resolution displayed on a heliocentric velocity scale. The profiles of Al III and Si IV have been scaled upwards by a factor of 70. The bottom 
panel shows the logarithm of the ratio of Al III and Si IV to S III as a function of velocity. The vertical dotted lines of the top panel represent 
the velocities of He? emission detected by Reynolds (1988b), though shifted by —2.2 km s _1 . 



